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ABSTRACT 


Flow control in packet-switched computer networks 
tries to avoid undesirable effects like congestion, 
deadlocks and also aids in improving performance in the 
dynamic sharing environment. Flow control functions are 
defined at various levels of protocol hierarchy* The 
efficiency of the scheme at a given level is limited 
by the efficiency of the schemes below it. At the lowest 
level are link control schemes, which have to account for 
random factors like noise on the links* Existing open 
queuing models are unable to account for the window effect 
while closed queuing models are unable to account for 
retransmissions due to nodal buffer overflow. A 2-stage 
queuing model, suitable for analysis of most schemes is 
presented in this thesis. The parameters of interest are 
calculated numerically. A simple approximation, to account 
for the underestimation of delay by Markovian assumptions 
under moderate traffic is introduced. Simulation results 
are also presented to test the accuracy of the model. 
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CHAPTER 1 


FLOW CONTROL IN PACKET-SWITCHED COMPUTER 
COMMUNICATIONS NETWORKS 

1.1 WHAT IS FLOW CONTROL 

In a dynamic sharing environment of packet switching 
close monitoring and control is necessary to ensure avoidance 
of undesirable effects like deadlocks and congestion. Routing 
[survey SCHW 80] and flow control [surveys KLEI 80 b,POUZ 8l] 
interact with each other m this task. At one extreme is 
loose coupling where flow control is activated only in case 
of a catastrophe. At the other end of the spectrum, is 
tight coupling where a packet is allowed to proceed only 
when the transit path is absolutely clear. 

Congestion avoidance requires two decision making 
processes : 

i) Choice of time-epochs at which a packet is admitted 
into the network (or a part of the network, e.g., 
subnetwork, a node etc.). The time epochs of admi- 
ssion at various levels of the network are controlled by 
mechanisms at various protocol hierarchies. For example, 
admission to network is controlled by Network Access 
Protocol while admission to next node, by link level 
protocol. 
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ii) Choice of path along which the packet reaches its 

destination with minimum delay subject to throughput 
constraints. 

It may be interesting to note that while (i) is 
sufficient for congestion avoidance, a combination of (i) and 
(ii) is necessary for optimal performance, i.e., minimal 
delay subject to throughput constraints. 

Routing accounts for decision process (ii). Decision 
process (i) is delegated to a set of mechanisms operated at 
various protocol hierarchies known as flow control. This 
also includes : 

i) error recovery mechanisms for data packets like re- 
transmissions or forward error correction (FEC). 
(Physical links for digital communication are better 
modelled by channels with memory (e.g. to account 
for error bursts) than the memoryless BS channel. 

A survey of the models, various error control schemes 
and their suitability is done by Karal and Sastry 
[KANA 78] ) # 

ii) Adaptive bit sampling rates or wordlengths for voice. 
This is discussed in Section 1.8. 

1.2 WHY PACKET SWITCHING ALLOWS BETTER CONTROL 

Availability of hardware at cheaper prices has made 
possible a number of alternatives to traditional circuit 
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switching* These include enhanced circuit switching and 
packet switching. An excellent comparison is done by 
Harrington [HARR 80], 

The nature of voice and interactive data traffic 
makes circuit switching inefficient. Voice traffic con- 
sists periods of talkspurt and silence. Interactive 
data consists of periods of activity and 'thinking time'. 

Enhanced circuit switching dynamically allocates and 
deallocates circuits on sensing activity or silence res- 
pectively. Time Assigned Speech Interpolation (TASI) 
[BUIL 59] was started on long trunk routes for voice and 
is found more efficient if the number of sharable channels 
is sufficiently large. 

The call connection is done by sensing beginning of 

a talkspurt. The fraction of traffic lost due to dynamic 

connection time or blocking due to all channels busy is 

known as cut-out fraction or freeze-out time. A cut-out 

fraction of 0.5 percent is acceptable for speech. Digital 

Speech Interpolation (DSI) is the digitized version. For 

interactive data, adaptive digital multiplexing, ADM, is 

used. This however does not perform as well as TAST 

because the error (cut-out fraction) tolerable is of the 
—7 

order of 10 for data. This factor also affects hybrid 
schemes as the fraction of data increases. 
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The only flow control measure in traditional CS 
systems is blocking. Voice traffic is better suited for 
this because of the delay constraints and no error recovery. 
Packetized voice is finding use. Weinstein [WE IN 78] showed 
that the cut out fraction in a TASI system is same as the 
fraction of packets dropped in a packet switched system with 
dropping if delay exceeds a certain threshold .Packet 
switching allows wider control than traditional circuit 
switching methods because 

i) Error recovery, advantages of adaptive routing are 
easily implemented in PS protocols. 

ii) Synchronization is crucial in CS systems. 

iii) Precedence and preemption are easily performed with 
PS protocols. In other words, PS systems can handle 
priorities better. 

iv) Packet switched systems can be designed for average 
load as opposed to peak load for CS systems. 

Enhanced circuit switched systems are able to incor- 
porate some of the advantages of packet switched systems. 

The distinction between PS and enhanced CS methods may be 
noted. If the number of outgoing channels is C, enhanced 
CS methods block subsequent customers if C are already 
active. PS methods allow more than C customers, with each 
customer facing greater delay. The mechanism of blocking 
can, of course, be introduced in a PS scheme as a call 
blocking protocol. 
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1.3 OBJECTIVES OF A FLOW CONTROL SCHEMA 

The main functions of flow control (FC) in packet 
networks are [KLEIN 80b]i 

i) prevention of throughput degradation and loss of 
efficiency due to overload 

ii) deadlock avoidance 

iii) fair allocation of resources amongst competing users. 

iv) speed matching between network and attached users. 

There are other objectives which are more related to a 
good economic balance between benefits and drawbacks of FC 
mechanisms [POUZ 81] , namely 

a) minimize response-time 

b) maximize useful throughput 

c) minimize transmission cost 

d) keep receiver busy at all times 

e) minimize overheads. 

The diversity of objectives indicates that there 
possibly cannot exist a single FC scheme. For example, 

( a) would be the most important factor for voice traffic 
while (b) and (c) may be more important for bulk data trans 
fer. As another example, objective (a) would require 
smaller packet lengths which would not satisfy objective(e) 

1 . 4 RESOURCES 


The resources required to carry a flow of data can be 
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categorized into two classes [POUZ 81] : 

a) basic j buffers, transmission bandwidth, processor time 

b) incidental : name space, table entries, logical channels, 

etc. 

Most studies consider buffers as the critical resource. In 
one case, CIGALE, mechanisms based on bandwidth were 
investigated. No study to the best of our knowledge has 
considered processor time a bottleneck. However, this 
aspect is expected to receive attention with faster 
communication lines e.g. T1 carrier ( 1.544 Mbits/sec) as 
compared to earlier 50 KBPS lines, coming into operation. 

Incidental throttling mechanisms are implementation 
dependent. There can be finite number of call establishments 
depending on the capacity allocated to the name space 
table. A packet may be blocked if its allocated logical 
channel (e.g. channel no = packet no. MOD no. of . channels) 
has not received the acknowledgement of previous packet 
even though the physical channel and receiver can accept 
packets. Also additional delay may be introduced by the 
overheads. 

1.5 PERFORMANCE MEASURES 

In a network there are two conflicting objectives : 

a) minimize delay 

b) maximize throughput 
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For example, throughput can be increased in a network 
by allowing more packets into it. But this means each 
packet encounters more waiting packets at each node. 

Hence, the delay is increased, and vice-versa. A performance 
measure that takes both effects into account is the through- 
put to delay ratio. This is termed ’power’. Bharath— Kumar 
and Jaffe [BHAR 8l] have done extensive work in the proper- 
ties of this function. One interesting aspect [JAFF 8l] is 
that flow control power is non-decentralizable. That is, 
any method to maximize power must consider the whole state 
of the network. 


This may be compared with the fact that there are 
efficient de-centralized routing algorithms. 

Lam [LAM 78] introduced a factor known as bursty factor. 
He argued that 'burstiness* of traffic on a network is not 
only dependent on the nature of generation at source but 
also on the characteristics of communication links being 
used. The factor is the transmission rate along the line. 


On a slow line, bursty data may be received more 
continuously whereas on a very fast line data generated 
with a high duty cycle can be received m a low duty cycle, 
and hence ’bursty*. 


Bursty factor 


fl ed message delay c o nstraint— 
Average "inter arrival time between messages 
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Bursty further is an upper bound of duty cycle. Traffic with 
low bursty factor is suited for PS. However, high bursty 
factor would not necessarily be better suited for CS. 

1.6 LEVELS OF FLOW CONTROL IN DATA NETWORKS 

Flow control can be exercised at various levels of 
protocol hierarchies [KLEIN 80b], The requirements are j 

a) Implementation of FC mechanism, at a level, must be 
consistent with other protocol functions at the same 
level. 

b) Interactions between various levels must be consistent. 

c ) There should be no duplication of functions at different 
levels. 

At the lowest level is the physical level. This is 
responsible for the electrical connections - activations/ 
deactivation. No flow control function is possible. 

Above the physical level is the link level. This 
layer purports to transport packets reliably across the 
physical links. Some of the flow control functions at this 
protocol level may be retransmission of packets not received 
correctly due to physical error, packets dropped due to 
congestion at next node, stop the sender (e.g. Receiver Not 
Ready (RNR) in HDLC, SDLC). HDLC (ISO defined standards) 

SDLC (IBM’s SNA), X.25 level 2 (subset of HDLC) are examples 
of link protocol implementations. 
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Above the link level we have the packet level protocol. 
This is used for establishing end-to-end user connections, 
through the network. This carries control information to 
route packets. The two methods of flow are : 

a) Virtual circuits (VC), and 

b) Datagrams. 

In datagram method, packets are controlled individually. 
Thus two packets frommthe same source may arrive at a 
destination node via different routes. In VC method, 
packets are received at the destination in the sequence 
of input. 

Datagram controls are based on : 

a) buffers, mostly , or 

b) bandwidth, e.g. CIGALE. 

In VC method the control is stream wise. Source-destina- 
tion pairs are identified. This method allows selective 
control on the various traffic classes. It allows one more 
degree of control by restricting the number of VC's. The 
controls can be : 

i) End-to-end - the inner sublayer (transit through 

intermediate nodes) is a datagram structure. Examples 
of these are in ARPANET, DATAPAC. 
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ii) Step-wise selective control is applied at node level 

as well. Examples are TYMNET, TRANSPAC, GMDNET. 

There are schemes where both datagrams and VC methods 
are used. These are known as hybrid schemes. An example 
is Integrated voice and data networks. Transmission of 
voice packets takes place in VC mode while data packets 
in datagram mode. 

Above the packet protocol is the end-to-end (ETE) 
protocol. The main objective of this is reliable trans- 
portation of single/multi packet messages from network 
entry to network exit nodes. Important functions of 
this protocol are reassembly of multipacket, messages at 
exit and regulation of input traffic using buffer allo- 
cation and windowing techniques. Some implementations 
do not have the ETE level in which case this function is 
delegated to higher level protocols. 

The highest level of network protocols which carries 
FC functions is the Transport level. This provides for 
reliable delivery of packets on the’ virtual circuits* 
between 2 remote processes. The goal is to regulate the 
flow so as to make the most efficient use of network 
resources. ‘Window* and Credit* schemes are generally 


used for this purpose 
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1.7 FLOW CONTROL METHODS USED AT VARIOUS LEVELS OF PROTOCOLS 
a ) Hop Level 

The main control techniques by observation of buffer 
sizes. Finer 1 s control may be exercised by observation 
on the number of unacknowledged packets (window size). 

The methods used are ; (i) Channel Queue Length (CQL) 

The classes of queues correspond to the various output 
channels. The restrictions placed would be the number of 
buffers each class can seize. At one extreme is the 
completely partioned (CP) scheme. This is dead lock-free 
but the utilization of resources is poor. The other 
extreme is completely shared (CS) scheme. The utilization 
is increased but this is prone to deadlocks and is unfair 
to input traffic of lower intensity. In practice a com- 
bination of the above methods SMXQ (maximum Q length), SMA 
(minimum allocation), SMXQA (maximum Q length with minimum 
allocation) are found to be good. Kamoun and Kleinrock 
[KAMO 80] have analyzed the performance of all the above 
s chemes . 

CQL avoids thruput degradation, unfairness, direct 
s tore-and-forward deadlocks. 

However, indirect S/F deadlocks cannot be avoided by 
CQL. This happens for example in a ring topology where each 
node wants to forward to the next node but each node has 


its buffer full 
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ii) Structured Buffer Pool (SBF) 

The buffers are divided into classes. The class of 
buffers that a packet can seize depends on the number of 
hops it has made. Indirect S/F deadlocks are prevented. 
This can be proved by resource graphs [KLEIN 80b]„GMD 
network uses this scheme [GIES 81]. 

iii) Virtual Circuits 

Selective control is made on each traffic class. This 
is the same as step-wise VC control scheme discussed in 
Section 1.6. 

b ) End-to-End 

This protocol is designed to work between source- 
destination pairs. The main objective is to protect the 
exit (receiver) node from congestion. An important by 
product is the prevention of global congestion, by reduction 
of the number of packets m transit in the network. This 
is based on end-to-end window scheme. 

c) Network Access 

The objective in this case is to reduce the transit de- 
lay at the cost of admission delays. A control on the 
total number of packets in transit m the network can 
reduce congestion. 
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The methods used are s 

i) Isarithmic Scheme [DAVI 79], Packets are allowed 
to proceed if they can secure a permit. The number of 
permits in the network is fixed. Isarithmic method 
requires an End-to-End scheme to maintain fairness in 
the distribution of the permits. 

ii) Input Buffer Limits 

Transit packets are given preference over entering 
packets. Also packets having travelled longer hops 
(closer to destination) are given preference over those 
which are fairly ‘newer* m transit. Lam and Reiser 
[LAM 79] and Kamoun [KAMO 8l] have studied these schemes. 

The number of buffers available to entering packets at 
gateway nodes, is restricted. 

iii) Choke Packet Scheme 

The network CIGALE uses this scheme. The emphasis 
is on the physical packet length rather than maximum 
packet length rather than maximum packet length in case of 
buffer schemes. This is done by observation of link 
loading. When the loading of a line sensed by a channel 
load limiter exceeds 0.7, a signalling packet, called choke 
packets is sent to sources generating excess traffic to 
thin link. If reduction of the loading does not take place, 
e.g. a source ignoring signal, a tripping mechanism auto- 
matically drops packets, on the link, causing congestion* 
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d) Transport Level 

The objective is to maintain maximum efficiency 
amongst processes. One of the factors to be considered is 
the loading of network to convey changes. For example, 
a graphics protocol (see discussion of 3 high level 
protocols, [SPR0 78]) is discussed by Sproullc and Cohen. 
Each change of scene would cause transmission of a large 
number of bits if information on all pixels are to be 
transmitted. By dividing the scene into segments, only 
the segment where change has occurred needs to be trans- 
mitted. 

1.8 FLOW CONTROL FOR VOICE AND INTEGRATED VOICE AND DATA 

NETWORKS 

The requirements for voice and data packets are 
somewhat complimentary. Data packets require high fidelity. 
This necessitates error recovery mechanisms. Delay is 
a secondary criteria to fidelity. Voice packets on the 
other hand, require a guaranteed delay, for continuity of 
conversations. Fidelity is not such a crucial criteria as 
loss upto 0.5 percent or slight distortions can be tolerated. 
The error recovery protocol is not activated for voice 
packets. Furthermore, to ensure guaranteed delay, voice 
packets are dropped if the waiting time at any node exceeds 
a threshold, chosen to satisfy the constraint of acceptable 
loss as well. Voice packets are generally sent m VC mode 
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while data packets can be sent in datagram or VC mode. 

Flow control of voice also involves applications of 
efficient speech coding techniques [GOLD 7 7]. Voice control 
methods have to be more adaptive. Voice controls are 
incorporated from end-to-end level of protocol hierarchy. 

As the network gets loaded, the bit sampling rate is 
reduced within acceptable limits. Bially et al [BIAL 80] 
presented an end-to-end adaptive control scheme based on 
embedded speech coding. Each frame of 20 ms is divided 
into pacicets of different priorities. Depending on the 
state of network, link management has to decide the cut off 
priorities. Packets below this priority level are dropped. 
This method can exercised to a point where the dropping 
ensures reasonable fidelity. This is advantageous over 
other schemes because of the ability of the link management 
to effictively reduce the bit rate. Other methods, involving 
end-to-end adaptive control suffer from the delay required 
to signal the end user. 

The delay restrictions and suppression of error recovery, 
made voice suited for circuit switching in hybrid PS and CS 
networks. Frame multiplexing structure [WEIN 80] is used. 

A frame of duration 10 ms is chosen m view of delay require- 
ments for voice. This is divided into two parts t 
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a) a fixed number of slots of fixed duration for voice 
transmitted in synchronous mode with circuit switching 

b) the remaining portion available for data m asynchronous 
mode. 

Fixed boundary case, when data packets cannot use vacant 
voice slots, is same as separate flow control schemes for 
voice, data traffic. When data is allowed to use vacant 
voice slots, the performance is better. This is the 
movable boundary case. 

TASI, DSI, packet voice control with dropping suffer 
from the problems of cut-out fractions. Adaptive bit 
sampling introduces distortions. Another control scheme 
for voice is Predicted Wordlength Assignment [GERH 82]. 

In this method the sampling rate remains constant. However, 
the number of bits per sample or wordlength, is adaptively 
varied. This method avoids cut-out fractions and distor- 
tions but introduces noise. 

1.9 FLOW CONTROL IN INTERCONNECTED NETWORKS 

Flow control in network interconnection, under virtual 
circuit mode of packet transmission (it may be pointed out 
the X.25 allows packet forward in VC mode) has been studied 
by Matsuto and Mori [MATS 81]. This is based on the virtual 
reduction of calls. 
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Edge [EDGE 79] studied the End point approach, which 
uses a standard internet end-to-end protocol and the hop- 
by-hop approach in which inter net calls consist of a 
sequence of local virtual calls across each intermediate net. 

Bennett [BENN 79] studied the problem of using trans- 
network architectures to support trans network bulk data 
transfer. 

1.10 FLOW CONTROL FOR REAL-TIME COMMUNICATIONS 

Cohen [COHE 80] pointed out the problems to be faced 
in flow control for real time communications. For example, 
in the schemes mentioned till now, a new packet is not 
accepted if the buffers are full. However, if the infor- 
mation carried by the packet is, say the position of an 
aircraft under surveillance, it would be useful to retain 
this packet and drop the older packet. The newer packet 
would convey more information to the receiver, in spite of 
dropping of intermediate information. 

He named this scheme 'Milk Rule' as opposed to the 
conventional ’Honey Rule’ where the latest packets to try 
entry are discarded. 

1.11 THE AIM OF THIS STUDY 

The efficiency of flow control scheme at a given level 
of protocol hierarchy depends on the schemes below it, at 
the lowest level is the link level control scheme. Problems 
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like noise in physical communication links affect the 
performance at this level. 

Voice control schemes generally exploit the inherent 
redundancy in speech signals by reducing the bit rate. The 
strmgest delay requirements achieved by dropping of 
packets at a node if a certain threshold is exceeded, 
further ensure that subsequent packets are not affected in 
terms of delay. 

Data packets , are more vulnerable to noise resulting 
in serious degradation of the network performance because 
of the additional traffic and delay introduced by the error 
recovery mechanisms. Improper choice of time out can cause 
redundant retransmissions. Improper choice of limit on 
unacknowledged packets in transit (window size) can lead to 
frequent blocking at next node and hence retransmissions, or 
lower thruput than what is possible. 

Consideration of these facts led to the development 
of a queuing model, the main result of this thesis, to 
study the performance of a general link-level protocol for 
continuous transmission with error recovery and finite 
window size. This model is found to be very useful for 
making design choice of the various parameters. 

Chapter 2 presents a summary of results and relevant 
applications that were motivating factors for the development 
of the model and the method of analysis. Chapter 3 presents 



19 


a step wise development of the model and a numerical pro- 
cedure to evaluate the necessary parameters. Chapter 4 
discusses the design of a simulation experiment, implemented 
using PASCAL on DEC-10, to study the accuracy of the model. 
Chapter 5 presents the results, observations and guidelines 
for further work. It also points out some modifications 
to the numerical procedure based on results derived in the 
Appendices, giving more accurate results for moderate 
traffic range. 
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CHAPTER 2 

PERFORMANCE EVALUATION IN COMPUTER NETWORKS 

2.1 INTRODUCTION 

Computer network nodes are modelled as queuing systems 
where the arrival process is the arrival of packets and the 
service process is the processing of the packets in buffer 
switching, transmission etc. These processes are strictly 
not deterministic. The packet arrivals do not occur at 
regular intervals of time and the packet lengths, which 
govern the transmission time (service) are not always fixed. 
Queuing effects become important in such systems where 
either the arrival process or the service process or both 
have variations. For example, when arrival and service 
processes are identical, and have variations, queuing 
theory predicts a non-zero expected queue length, contrary 
to intuitive thinking. Results of queuing theory are used 
to predict the relevant measures like average delay at 
nodes, throughput, buffer overflow probability, etc. 

Queuing theory was used by Erlang to study the per- 
formance of telephone systems. His famous Erlang B-formula, 
the probability of C customers in M/M/C/C queuing system, 
gives the call blocking probability and is used extensively 
in telephony. 
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Kendall [KEND 51] pointed out the three basic factors 
in queuing systems : 

a) arrival process 

b) service process 

c) service discipline. 

The shorthand for representation of queuing systems 
was introduced by Kendall [KEND 53] • 

There are two approaches to the analysis of queues ; 

i) difference-differential equations of Erlang 

ii) integral equations of Lindley [LIND 52] . 

2 . 2 SIMPLE QUEUES 

The simplest queuing system is the M/M/l/oo system. 

The arrival process is poisson and the service time 
density function is exponential. The poisson arrival 
process gives rise to exponential inter-arrival time 
density function. Observations show that the demands of a 
large independent population, e.g. telephone call demands, 
tends to follow poisson distribution. 

The exponential density function is the only one that 
gives rise to memoryless condition, 

P[u < t < x + 6] = P[t< 6] 

-Xt 

Let the inter-arrival time density function be, a( t)= X e “ ,t>o. 
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and the service time density function be, 
b(t) = p.e"’^ it , t > 0 

Then, the state of the system is completely specified by 
the number of customers in the system whose probability 
[ KLEI 75], 

p n = ^(l-P), n = 0,1,2, (2.2.1) 

where P = X /p. 

The condition for stability is P < 1* 

Computer Network nodes have finite buffers. If N is 
the maximum number of customers in service or waiting, 
then the system M/M/l/N gives rise to the state probabili- 
ties [KLEI 75, SCHW 77] 

r (i-p)p n /(i-p N+1 ), p ,4 i 

p n = J (2 - 2 - 2) 

l/(N+l) , P = 1 

The ergodicity condition is P < °® . 

The probability that an arriving customer is turned 
away due to no waiting room is the probability that N 
customers are already waiting or in service at the system. 
This is given by 


= (l-p)p N /(l-P N+1 ) 


(2.2.3) 
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The expected number in the system, L, is given by [GROS 74] 


i) 

ii) 


L = 
L a 


P/( 1-p) for M/M/l/o» 
f p [ i~(n+i)p n + np n+1 3 
j (i-p)d-p k+i ) ’ 


system 
P £ 1 


^ N/2 , p = 1. 

for M/M/l/N system. 


(2.2.4) 


The waiting time in the system (including service) 
is given by Little’s Formula [LITT 61, JEWE 67] : 


i) W = L/X , 

ii) W = L/ x * , M/M/l/N (2.2.5) 
where x ' = X (1-p^) (2.2.6) 
is the effective throughput. 


Apart from finite capacity queues, variations of the 
M/M/l system applicable to computer networks include : 

i) State dependent inter-arrival times and/or service 

times. For example, the probability of arrival of a 
talkspurt or length of a talkspurt (arrival of 
silence period) depends on the number of active 
talkers [WEIN 78], The analysis of W-R retrial model 
for telephone systems [FRED 79] incorporates the 
effect of a constant retrial probability by blocked 
customers by a state dependent service rate, by taking 
the marginal probabilities. 



24 


ii) Queues with impatience applicable, e.g., to voice 
packets which are dropped after a certain delay 
threshold at a node [WE IN 78] . 

iii) Bulk arrival or service, e.g. the arrival of acknowledge- 
ments under IBM’s SNA virtual pacing scheme can be 
modelled as a bulk arrival process [SCHW 82] • 

iv) Queues with finite population. Physical computer net- 
works or systems cater to a finite population. However, 
if this is sufficiently large the finite population 
constraint, modelled via state dependent methods, is 
dropped. An application of finite population model 

is finding the number of interactive terminals attachable 
to a system beyond which the response is greatly deterio- 
rated [KLEI 76]. This is known as the saturation number, 
M*, given by, M* = 1 + Jj- • 

Analysis of these cases of the M/M/l model is available 
in most books on queuing systems, e.g. [KLEI 75]. 

2.3 QUEUES WITH INTERMEDIATE COMPLEXITY 

The queuing system becomes more complex when one of the 
processes is not Markovian. These are the M/G/l queue where 
the service time can take a general distribution and Gl/M/l 
queue where the inter-arrival time distribution can be 
arbitrary. 
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2.3.1 M/G/l Queue 

The analysis of the M/G/l queue uses the embedded 
Markov process approach of Kendall [KEND 51, KLEI 75]. The 
state of the system must now include the service already 
expended on the customer in service. For the M/M/l queue 
the incremental service time required for completion for 
the customer in service is independent of the service 
received due to the memoryless property. The state 
vector [N(t), X(t)] where N(t) = number in system, 

X(t) = service time spent on customer in service, is now a 
Markovian process. Kendall’s ingenuity [KEND 53] was to 
consider the embedded Markov process, the state of the 
system at departure epochs. 

The state equations can now be written as, 

3+1 

Pj = I P]c d j+l-k (2. 3. 1.1) 

k=0 


where d^ = probability of k arrivals since last departure. 
If b(t) is the density function of the service time and , 
the parameter of the poisson arrival process, then, 


d 


k 



Xt (X t) k 

lei 


b(t) dt 


(2. 3.1.2) 


The parameters of interest can be found by transform methods 
[KLEI 75], The z-transform of the system queue size is 
given by, 
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(1- p)(l-z) 

Q(z) = B*( X - Xz) (2. 3. 1.3) 

B *(X -X z)-z 

where B*( s ) = Laplace Transform of service time pdf and 
is the parameter of the poisson arrival process. This is 
known as the Pollaczek-Khintchine transform equation. Taha 
[TAHA 71] has given the results for M/E^/l and M/D/l systems. 

The Laplace transform of the waiting time in queue is 
given by 

W*(s) = Sll=J L l (2.3.1.4a) 

s- X+ B*( s) 


The overall system waiting time is given by S*( s ) = W*(s)B*(s) 

(2.3.1. 4b ) 


Finch [FINC 58] analyzed the M/G/l system with finite waiting 
room of size, N, He showed that if ( l-x)/(d( x)-x) can be 
expanded as a power series in some region, then the generating 
function for state probabilities of interest is given by 


P n ( 1-x) 

P(x) =-2 (2. 3.1. 5) 

d(x) - x 

where P(x) = P Q +P 1 +P 2 x + ... + P^" 1 + ? N+1 X N , 


P n = probability that number in system at departure 
epoch = n, 


and 

d(x) = l d k X k 


k=0 
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where is given by eqn. (2. 3. 1*2), P Q is determined to 
ensure that the sum of the probabilities is equal to 1. 

Courtois [COUR 80] analyzed the M/G/l finite queuing 
system in a different manner. The open M/G/l/N queue is 
replaced by a closed two stage cyclic queuing network with 
N cycling customers with exponential server (arrival process) 
at first stage followed by a general server at the second 
s tage. 

Considerable work in computer systems, e.g. scheduling, 
has been done under M/G/l consideration. Multiprogramming, 
time sharing systems are advantageous when the coefficient 
of variation of service time is greater than 1 [KLEI 76]. 

The analogy to networks is that if the coefficient of 
variation of message lengths is greater than 1 packetization 
gives better performance. 

A stop-and-wait protocol for link level was analyzed by 
Fayolle et al [FAYO 78] under M/G/l assumption. Their obje- 
ctive was to predict optimal time-out periods to s 

i ) minimize delay 

ii) maximize throughput 

iii) minimize buffer overflow probability with finite 
buffers • 

Kuehn [KUEH 81] analyzed the performance of two 
different ARQ-protocols for half duplex transmission between 
a primary station and an arbitrary number of polled secondary 
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stations. He analyzed them under cyclic queuing models 
of the type M^ x Vg/ 1 with poisson batch arrivals and general 
service time processes. 

2.3.2 GI/M/1 Queue 

Let the arrival process have density function a(t) 
and service time density function be exponential with 
parameter p. 

The observation epochs are now the arrival instants. 

We have now, the state probabilities, 

p j = lo b+lc-l d k (.2.3.2.!) 

where d^, the probability of k customers services since last 
arrival, is given by, 

« -ut , ,sk 

d. = / L&kL. a(t) dt (2. 3. 2. 2) 

Ki 

The infinite summation in equation (2. 3. 2.1) as opposed 
to finite summation for the M/G/l case (eqn. (2. 3. 1.1)) leads 
to a simple results for the state probabilities [TAHA 71, 

KLEI 75], 

p. = (1 -cr)a j (2. 3,2.3) 

where, a is the solution in (0,1) of 

cr = A*(p-pa) (2. 3.2. 4) 

where A*(s) is the Laplace Transform of a(t). 
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The waiting time for Gl/M/l queue follows exponential 
distribution. 

An exponential server fed by the output of an M/G/l 
system constitutes a Gl/M/l system,. 

2.3.3 The M/G/l Conservation Law and Fairness Policies in 
Scheduling in Computer Networks 


In computer networks, priorities have to be assigned 
to traffic classes. It is useful to find a way for evalua- 
tion of such systems. Kleinrock’s famous M/G/l conserva- 
tion law states [KLEI 76] that for any non-preemptive work 
conserving queuing discipline. 


y p w = 

4 M p p 
p=l 


PW 


o 


1- p 


, p <1 

» P >1 


( 2. 3* 3 • l) 


where is the traffic intensity of the pth class and 


p - p • W is the mean response time for the pth class 


P-1 


and W Q = mean response time (service + waiting) using the 
Pollaczek-Khintchine formula. 


The conservation theorem allows the network designer 
to reallocated the delays among classes without significantly 
affecting the overall mean-end-to-end delay. 

Wong et al [WONG 82] applied this to the study of 
fairness. They defined a generalized fairness measure, 
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P f y 

F = -py* l ~~r (Tp - Tp) 2 (2. 3,3. 2) 

T P=1 

where, T = mean overall delay 

Tp = mean delay for pth class 
+■ 

Tp = target delay for pth class 

Yp»y are throughput for pth class, overall traffic 
respectively. 

Let = a K (2, 3. 3.3) 

Jr r 

where a is a parameter indicating the priority of the pth 

Jr 

class and K, a constant that can be found by application of 
the conservation law, 


P 

l 


P =i 



Therefore, 


K - 


1 T 

p 

pii 


Y CL 
J p p 


(2. 3. 3.4) 


(2.3. 3. 5) 


The various priority schemes considered were : 

i) uniform delay fairness s a = 1 for all p. 

ii) Hop dependent fairness ; This case is applicable for 
location independent tariffs (e.g. Telenet). One might 
want the delay of a class to be proportional to the 
number of hops on the path used by that class. 
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This is done by setting 

a = 7c } the number of hops. 

P P 

iii) Traffic dependent fairness : Higher priority for class 
with higher traffic volume is allocated. 



FCFS scheduling policy gives approximately the same 
fairness measure as the optimal policy for the hep dependent 
case. For the other two cases, FCFS performs poorly. 

2.4 COMPLEX QUEUING SYSTEMS 

The most complex queuing system is Gl/G/l system where 
the arrivals and service times can take arbitrary distri- 
butions. Use of Gl/G/l methods would allow prediction of 
upper bounds on delays in computer networks. 

The method of analysis for Gl/G/l queue uses a random 
variable. 


u_ = x „ 
n n 


n+1 


(2.4.1) 


where x„ = service time for nth customer . and 
n 

t n+ ^ = inter-arrival time for the (n+l)th customer. 


If A*(s) and B*(s) are the Laplace Transforms of inter- 
arrival time and service time density functions, then the 
Laplace Transform of the density function of the random 
variable, u, is given by [KLEI 76] 
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C(s) = A*(-s) B*(s) (2.4.2) 

The distribution of the waiting time is given by Lindley* s 
Integral equation [LIND 52, KLEI 76], 

r Y 

w(y) a J / W(y-u) dC(u), y 0 (2.4.3) 

J 

1^0 , y < 0 

Kingman [KING 62, KLEI 76] showed that the limiting waiting 
time distribution under heavy traffic (traffic intensity 
close to 1 but strictly less than 1 to maintain stability) 
is exponential with mean 

(0 a + °b )/2( e) * (2.4.4) 

where , 

a a , a b - standard deviation of inter-arrival, service 
times respectively 
p = traffic intensity 
t = mean inter-arrival time. 

Recent developments in Gl/G/l queues include : 

1) analysis by discrete fourier transform by Ackroyd 
[ACKR 80] 

2) approximation of the waiting time distribution, W(x) in 
the mean and probability of delay, P d = 1-W(0), by 
exponential distribution of the form, 

W Jx) = 1-Ce~ ax , by Fredericks [FRED 82]. 

3 . 
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2.5 APPROXIMATION METHODS 

Fluid approximation (approximation in mean) and diffu- 
sion approximation (approximation in mean and variance) are 
useful in performance evaluation [KLEI 76, KOBA 74a, b] 
especially since the exact methods are difficult for com- 
plex systems. These methods are applicable under moderate 
to heavy traffic. This is the range where the effectiveness 
of flow control schemes lie. 


The above methods replace the discrete stochastic 
process by a continuous process. The basis is as follows : 

If a(t) = number of arrivals In (0,t) , and 

6(t) = number of departures in (0,t), 

then, if a(t) is largo compared to unity, 



(2.5.1) 


The same holds for 5(t), 6( t ) . Thus, the discrete 
stochastic processes a(t) and 6(t) are replaced by the 
continuous deterministic processes lx( tj and &( t) for fluid 
approximation. The number in the system, 


N( t) = aTTI - 6( t) (2.5.2) 

The diffusion approximation method replaces the dis- 
crete stochastic processes a(t), 6(t) by continuous 
Gaussian distributions by accounting for the variations as 
well [KLEI 76]. 
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2.6 OUTPUT PROCESSES AND NETWORK OF QUEUES 

Output processes are important because not all system 
parameters are observable and a measure of performance is 
throughput [PACK 78], Output processes are also important 
when extension to network of queues is considered. Burke 
[BURK 56] showed that the output of a M/M/l system is 
poisson. Finch [FINC 59] pointed out that for the M/G/l 
system the interdeparture times are independent only if 
the distribution of service time is exponential. For 
finite waiting room cases, even with exponential server the 
independence is not maintained. Daley [DALE 68] showed 
that the outputs of Gl/M/l queue are independent if 
arrival process is poisson. Mirasol [MIRA 63] showed 
that the output of a M/G/» queue is poisson. The 
engineering implication is that if the waiting time at a 
node is sufficiently small (effect of a server always 
available), poisson output may be assumed. Recent work 
on output processes include a survey by Burke [BURK 72] and 
some results on M/D/s queue [PACK 78]. Burke [BURK 72] 
pointed out that the output of the M/M/l/N queue, when the 
dropped customers are also included, is poisson. 

These factors restrict closed form solutions to 
network of queues to limited cases. 
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Jackson [JACK 57] solved the problem for a network of 

M/M/c/*® queues. He arrived at the famous result known as 

Jackson’s theorem. The joint probability of n^ customers 

at service station 1, . n^ customers at service station 

i etc. is the same as the product of the probabilities of the 

state at each node, independently. Or, 

M 

p(n 1 ,n 2 , n m ) = R Pj(n.) (2.6.1) 

3=1 

Burke [BURK 76] proved a fact conjectured earlier. The 
output of a M/M/l node with feedback is not poisson, in 
spite of the fact that the queue size, at each node of the 
Jackson network, is as if it were fed by independent 
poisson input. Kobayashi and Reiser [KOBA 75] pointed out 
that the factor of interest is actually the mean work demand 
brought by a job, and not how this work is divided into 
individual visits to a node by routing. 

In computer networks the assumptions of Jackson's 
theorem are not satisfied because [KLEI 76] : 

i) service times are proportional to message lengths and 
can be anything but independent 

ii) arrival times are conditioned on service times at nodes 
and thus cannot be independent. 

If, however, the number of input channels to the 
node as well as the number of outgoing channels are sufficiently 
large then the dependencies reduce* It is also known that 
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- N_ - V 

requests from a set of large number of independent sources 
to a pool of shared resources form a poisson process. 

These factors led Kleinrock to formulate the famous 
Independence Assumption for Computer Networks, The basis 
for service times at different nodes being exponentially 
distributed independent random variables was thoroughly 
tested out by simulation. 

In a packet network of fixed size packets, it would 
appear that service times would be constant. However, 
considering overall service factors like processing time at 
next node, loss, delay in receipt of acknowledgement etc,, 
the exponential distribution assumption is a fairly good 
one, 

Gordon and Newell [GORD 67a] showed that under certain 
conditions closed queuing networks lead to product form 
solutions. They also considered the finite queue case 
[GORD 68b], Using the results of [GORD 67a], Buzen 
[BUZE 73] presented algorithms for efficient computation, 
Kuehn [KUEH 79] used approximations by decomposition method 
to analyze network of queues under Gl/G/l assumptions. 
Kaufman et al [KAUF 81] used sparse matrix methods for 
analysis using Markov state space model. More details can 
be incorporated, e.g., deadlock probability in suppressing 
acknowledgement. However, the complexity of solution in- 
creases as the number of states increase. 
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The g oncralization of Jackson’s theorem came from 
Baskett ot al. [BASK 75]. They showed that under the 
following conditions : 

I ) server exponential - discipline i’CFS. Service time 
distribution is identical for all classes of customers. 

II ) Processor sharing, server having rational Laplace 
Transform. Each class m~>y have distinct service time. 

iii) The number of servers is greater than the maximum 
number of customers. Each cl^ss of customer may have 
distinct service time, having rational Laplace Transform. 

iv) Discipline LIFO, servers having rational Laplace 
Transform. Each class may have distinct service time. 

closed form solutions result. These systems are known 
as BCtoP queues, 

Cox [COX 55] showed that any distribution having rational 
Laplace Transform can bo represented by stages of exponential 
servers has shown in Fig. 2.6.1. This includes feedback systems. 
Another representation can be given m terms of parallel 
exponenticl stages. 

The conventional approach to the solution of queuing 
networks makes use of global balance equations, Whittle 
[wHIT 68] introduced a simplification known as independent 
balance equations. These equate the rate of flow into a 
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state by a customer entering a stage of service to the 
rate of flow out of that state due to a customer leaving 
that state of service. Each global balance equation is the 
sum of independent balance equations. 

Discrete time queuing systems and their networks 
have been studied [BHAR 80, HSU 76], Examples of tnese 
include modelling of processor time (service) which is a 
multiple of a fixed quantum or allocation of slots for 
transmission which are multiples of a fixed quantum. Inter- 
arrival times may also be discrete, in terms of multiples 
of a fixed quantum. Unlike the Jackson’s analog network, 
a network in which packets arriving in independent Bernoulli 
stream, are being serviced by a resource with geometric 
distribution, the independence property of outputs is not 
maintained. 
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2.7 APPLICATIONS OF NETWORK OF QUEUES 
2.7*1 Some Open Queuing Models 

Schweitzer and Lam [SCHWE 76] used the results of BCMP 
queues to calculate the blocking probability at a node with 
many outgoing channels under sharable buffer scheme, under 
markovian assumptions for arrival and service times. 

They gave a closed form solution to their model of 
the form, P (n Q , N, M, L ) , where n Q = number at CPU (input) 
queue, N = (n^, ..., n m ) , where n^ is the number at service 
center i for transmissions. M = (m^, m^) where irn is 

the number transmitted and waiting for acknowledgement, and 
will be acknowledge at service station i. L = (1^, l m ) 

where 1^ is the number transmitted but will not be acknowledge 
and will be fed back to end of transmission queue after time 
out. 

It is assumed that once a packet is transmitted there 
is always a server available to handle the acknowledgement 
process (This allows application of assumption 3 of BCMP 
queues ) . 

Irland et al [IRLA 77] gave a numerical method of 
analysis for a network of nodes in tandem. The influence 
of finite buffering at each node, on blocking input paakets 
and their subsequent retransmissions was considered. 
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Their results show that host-retransmission schemes 
and node retransmission schemes result m approximately the 
same delay, under the assumptions made. 

2.7.2 Some Closed Queuing Network Models 

Window schemes are easily modelled in terms of closed 
queuing networks. A finite number of customers (window 
size), can be circulating. This is because at one extreme 
W packets can be traversing and at the other extreme 
the W acknowledgements can be traversing back, and 

a combination summing to W where W is the windows ize, is 
possible. 

Reiser [REIS 79] modelled a network using window 
control. The open queuing network was converted into a 
closed multichain queuing network. For example, chain 
could be due to end-to-end window. Another chain could be 
due to the link-level window. Mean value analysis of 
Reiser and Lavenburg [REIS 80] and convolution method 
[REIS 81a] are used to analyze these networks. 

Schwartz [SCHW 82] showed that the optimal end-to-end 
window size for a L-hop path Is, W = L. This assumption 
assumed identical servers. This allowed replacement by a 
single state dependent server by use of Norton’s theorem 
[CHAN 75]. Also, the analysis was performed under satura- 
tion conditions (traffic intensity = 1) on the basis that 
this is the region where effect of flow control is most 
prominent. 
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Reiser [REIS 8 lb] introduced the effect of admission 
delays in a VC route caused by blocking due to lack of 
permit (window full) and showed that the optimal window for 
the end-to-end scheme is, now, 

k 0 = 2L 

End-to-end data flow analysis based on closed queuing 
networks with application to HMINET was done by Baum et al 
[BAUM 79]. They found that the size of buffer space needed 
for transaction buffering and throughput are dependent on 
workload of the local comouter. They thus, recommended 
a variation of the level of multiprogramming as a possible 
traffic control measure* 

2.8 ZERO INSERTION AND ITS INFLUENCE ON PACKET LENGTHS FOR 

HDLC PROTOCOLS 

The general format of an HDLC information (I) frame is 
shown in Fig. 2.8.1. 


F 

A 

C 

T 

I 

i : 

. PCS 

hi 


Fig, 2,8.1 


It consists of 

i) a one byte flag ( F) , 

ii) a one byte address field (A), 

iii) a one byte control field (C) 
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iv) arbitrary bytes of information (I), 

v) a 2 byte frame check sequence (FCS),and 

vi) the one byte flag (F), 

The other formats e.g. supervisor are of fixed length. 

It is the I frames which have major influence in the per- 
formance. 

The one byte flag (F) consists of a single *0* followed 
by 6 ones and terminated by a ’O’ - 01111110. This is used 
for synchronization. To avoid false synchronization due to 
this sequence in other fields, a 0 is added after each 
sequence of 5 consecutive ones. Ma [MA 82] derived the 
results for efficiency under general conditions where, a zero 
is added after n consecutive one’s, using a Markov model. 

The optimal value of n, is, 

n opt = log 2 { 2 + + ^ (2+ - 4} -1 (2.8.1) 

where L is the frame length. 

The efficiency for general n and length L, 

EFF(n ’ L) ° F T CnTL'I (2 * 8 ’ 2) 

0 n+l 1 

FL = L • ~ r= ^ + 2(n+3) (2.8.3) 

2 n+1 -2 

where FL = resultant length due to zero insertions. 
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Ma’s results show s 


i) under noiseless lines EFF(5,L) is within 1-2 percent 

of EFF L) hence the HDLC choice is good 

ii) under noisy lines, EFF ( 5,L) is not as good as 

EFF ( n 0 p£»k)» beyond a limit for L. This is because 
greater the frame (packet) length the larger the 
number of zero insertions. One solution is to restrict 
the packet length to the acceptable value of L, and 
use the present HDLC format. Ma, on the other hand 
suggests a flag field of 2 bytes and FCS of 4 bytes. 
This is to make the flag field consisting of a ’0*, 

14 consecutive I’s a 'O*. Now the zero insertion takes 
place after every 13 consecutive l's in other fields. 
This is more efficient and EFF ( 13, L ) is found to be 
closer to EFF ( n 0 p^>^) for L exceeding certain 
limit. 


2.9 DYNAMIC FLOW CONTROL 

Kermani and Kleinrock [KERM 80] studied dynamic flow 
control. The parameters like time out, window size under 
end-to-end control are chosen by a Markovian decision model 
based on previous observations. This is found to give better 
performance than static flow control [KLEI 80a]. 

2.10 FLOW CONTROL IN INTEGRATED VOICE AND DATA NETWORKS 
Analysis of flow control in integrated voice, data 
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networks are also being carried out. Fischer and Harris 
[FISC 76] presented a model for analysis of the integrated 
multiplex structure where each frame in subdivided structure 
where each frame m subdivided into fixed slots for voice in 
CS mode and the remaining for data in PS mode. 

Weinstein et al [UEIN 80] studied by simulation flow 
control of voice and data. Data packet control is effected 
by buffer limitation while voice packets are controlled by 
adaptive bit sampling. Smoothing algontnms [JAYA 76] can be 
used to handle transmission errors for voice. 

Ross and Mowaffi [ROSS 82] presented an M/G/l model 
for study of the parameters of data packets in a hybrid 
scheme. The server is modelled as one with rest periods 
(to service the voice slots in CS mode). The rest period 
is constant for a fixed boundary case. 

2.11 PROPOSED QUEUING MODEL 

Open queuing models are useful for predicting throughput, 
delay and buffer overflow probabilities, but do not account 
for blocking due to windowing effect. Closed queuing models, 
on the other hand, do handle window effect. However, the 
buffer overflow probability and retransmissions are not 
handled by these models. 

This led to the development of a model for link level 
flow control incorporating buffer overflow, retransmission 
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after time out and windowing, presented in Chapter 3. Open 
queuing models of a node are useful for extension to network 
of queues and finding suitable approximation using BCMP 
results. Using the ideas of Gordon and Hewell [GORD 67a], 
transformation of a closed queuing system into open queuing 
system with state dependent input buffer allocation is 
done, A numerical method, using ideas of Irland et al 
[IRLA 77] and Hillier et al [HILL 6 7], with considerable simpli- 
fications is formulated. It is found that for most acceptable 
values of buffer size and window in networks the input buffer 
allocation can be assumed independent and hence the open 
model is found quite accurate. 

An easy to evaluate approximation is introduced to 
account for the under-estimation of delay with M/M/l assumption 
under moderate traffic. Use of the exact methods, e.g. M/G/l/K 
would be far more complex. 
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CHAPTER 3 

A QUEUING MODEL FOR SOME FLOW CONTROL SCHEMES 

3.1 INTRODUCTION: 

In this chapter, we develop a queuing model for link - 
level flow control schemes incorporating retransmissions, 
time out and window. This model is depicted in Fig. 3.1. 


Time out 



Assumptions made m the model are: 

Arrivals to the node are poisson with parameter \ . The 
packet lengths are distributed according to exponential 
distribution. The time for transmission proportional to packet 
length is exponentially distributed with parameter This 

is represented by the server at stage 1, 

At most N packets can be presented at the node. This 
corresponds to physical buffer of size, N. A transmitted 
packet continues to occupy the buffer till receipt of 
acknowledgement. 
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The buffer management scheme tags this buffer to the 
list of unacknowledged packets. Each packet is acknowledged. 

The acknowledgements are m order of transmission of packets. 
The process of waiting for acknowledgement after transmission 
is represented by stage 2. The queue, here, represents the 
fact that the acknowledgement for a packet is sent after the 
acknowledgement for the previous packets. The size of this 
queue (maximum number waiting or in service) is Vi, the 
window size. 

The acknowledgement intervals (time between two acknow- 
ledgements) are assumed exponentially distributed, with parameter 
^ 2 * to account for the various factors like delay m processing 
of packet at next node for detection of transmission errors 
and despatching acknowledgement, piggybacking of acknowledgements 
on packets m reverse direction, propagation delays etc. It 
may be noted that the transmission time for a packet is the 
elapsed time between sending the first and the last bit of 
the packet. 

The buffer management scheme explained earlier, leads 
to an arriving packet at time t, to the physical node finding 
the buffer availability for packets at the transmission process 
(stage l). 


N* (t) = N - w( t) 


(3,1,1) 
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where w( t) is the occupany by packets waiting acknowledgement 
at time t. Thus, the inherently closec' queuing system with 
i J circulating customers^ converted to an open queuing system 
with the buffer availability N*( t) dependent on occupancy 
by unacknowledged packets w(t)» The analysis of this system 
would be complex, we shall simplify further, shortly. 

If all N*( t) buffers are full (corresponding to all N 
buffers and full) an incoming packet is dropped. Let the 
probability of this event be p^. 

Consideration of physical buffer sizes, window and 
expected occupancy by unacknowledged packets led us to make 
a simplification. The average buffer occupancy by unacknow- 
ledged packets is small as compared to N. Thus, we have the 
approximation, from (Eqn. (3.1.1)), 

N*( t) = N ( 3.1.2) 

We have, therefore, modelled a completely shared buffer 
management scheme as a completely partioned scheme with N 
buffers available to incoming packets and Y.< for unacknowledged 
packets. Further more the complex system is now modelled as a 
tandem 2-stage queuing system, with independent, finite buffers. 
The average accupancy by unacknowledged packets is of the order 
of 2 to 3 (this is confirmed by results of simulation), when 
the acknowledgements arrive at the same rate as transmission of 
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packets. This is because the finite buffering (stage l) 
blocks input traffic to physical node under heavy traffic 
conditions. The acknowledgement process, therefore, caters 
to an already regulated traffic. The acknowledgement process 
is much faster than transmission process under normal 
conditions. This is because the acknowledgement packets 
are much smaller, and due to the effects of acknowledging 
upto a certain number of packets in a single acknowledgement, 
piggybacking etc. 

The above arguments would require careful review for 
satellite and long-distance terrestial channels. The window 
size has to be much larger (e.g. m HDLC extended mode, 

V'J = 128) to allow continuous transmissions m presence of 
larger propagation delays (round-trip delays). 

An accepted packet receives service at stage 1, which 
corresponds to the transmission process. The queue discipline 
is FCFS. With a probability, p e , the packet has to be retrans- 
mitted. This may be because of erroneous delivery of packet 
or acknowledgement m transmission or blocking due to buffers 
full at next node. The retransmission takes place after a 
fixed duration, called time-out. This procedure repeats till 

successful transmission. iNe assume that time-out = l. mean 

-■1 

transmission time = i . 77 ™ (3.1.3) 

F-i 
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The Go-Back~N method of retransmission is used because 

of the sequencing assumptions, A packet will not be accepted 

out of sequence i at a node. In the actual flow control 

scheme, when a packet has to be retransmitted after the timer 

interrupt (time-out), the subsequent packets (at most W-l) 

are also retransmitted. If i> W then the server at stage 1 

is idle on the average for the duration, 

Time out - time to transmit ''! packets, 

if an unsuccessful transmission takes place. In our model, 

we incorporate this effect by assuming that server 1 is 

busy during the entire time-out period and include this time 

in the service time of the packet requiring retransmission. 

The transmissions of other packets in-between are ignored. The 

constant time-out period is modelled as an Erlang-i stage 

•» 

with i exponential servers m series, with mean , each. 

With the above considerations the average service time 
of a packet is. 



H riMsaur — j-«-— timsovt — ^ 


X Qff = (l-p e )jC+p G (l-p e )(i+2)X+ P g (l-p G )(2i+3)X+ 


+ p e J ( l-p e )[j( i+l)+l]X + 


(3.1.4) 
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where X = 1 f\x^ = mean transmission time. 

The transmissions of packet number 2 due to retransmissions of 
packet 1, have been futile. Effectively packet 2 is treated 
as starting a fresh attempt to transmit (after 1 ) but with 
corresponding increase m its waiting time. After successful 
transmission, the packet moves over to the queue of unacknow- 
ledged packets m stage 2 of the model. 


In the flow control scheme using a window of size ■/, 
transmission of fresh packet stops after h unacknowledged 
packets have accumulated. In our model this is achieved by 
assuming that only # buffers are available at stage 2. Server 
1 gets blocked when all these buffers are occupied. Let the 
probability of this event be p 2 « 


The delay faced by a packet m the system is the delay 
from arrival till departure from service m stage 1. This 
represents the delay transmission completion of the packet 
to Hie next physical node. The second stage represents the 
events till the copy of the packet, kept for retransmissions 
if necessary, is discarded. The role of the second stage is 
in blocking the first stage and hence changing its rate. The 
extent of blocking depends on the flow control scheme e.g. 
for stop-and-go, P 2 =l ? while for continuous transmissions (no 
window constraints) pg = 0. 


CENTRAL JJSRMCT 

rCTsrm.- 
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The model as formulated is quite complex and its 
analysis is not straight forward. To simplify the task we 
carry out the analysis m a step-by-step manner. We do this 
by making various simplifying assumptions which reduce the 
model to easily-solvable cases. These cases correspond to 
simpler versions of the generalized flow control scheme. 
Their analysis helps us in finally analyzing the general 
model. 

The various versions considered by us ares 

a) Continuous transmission with infinite window 

1} no error recovery 
li) retransmission after zero time out 
111) retransmission after finite time out 

b) Send-and-wait schemes 

c) Continuous transmission with finite window, no 
error recovery and, finally 

d) Continuous transmissions with finite window and 
retransmission after finite time-out (most 
general version). 

The remaining sections of this chapter describe the 
analysis of these versions. 
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3.2 CONTINUOUS TRANSMISSION *'ITH INFINITE WINDOW; 

Assuming an infinite window size implies that the 
probability P 2 of blocking due to second stage is zero. Since 

the second stage influences the first stage only through block 
blocking , we can eliminate the second stage. The system is 
now a single server with finite waiting room which is easy 
to analyze. 

3.2.1 No Error Recovery; 

This is equivalent to assuming that p e =0, so there are 
no retransmissions. (There may be retransmissions due to 
acknowledgements delayed beyond time out. Analysing the 
effect of these redundant retransmissions is postponed to 
avoid complication). Hence the service time distribution 
is simply exponential with parameter, [i^j reducing the 
entire system to M/jvl/l/N whose performance can easily be 
determined (Sec. 2.2). 

3.2.2 Retransmission after Zero Timeout; 

This corresponds to a hypothetical system m which 
the acknowledgement information is available as soon as 
it is transmitted. Note that the analysis of servers m 
series is done m time domain by convolution, which corres- 
ponds to product in L-T domain. 
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Referring to Fig. 3.1, with Time-out = 0, the Laplace 
Transform of the service time, 


Pi 


4 


1 \2 


BfCs) - ( 1—p e ) ( s _|_| i ) + P e ^ l“P e ^“Pe ^s+p^ 


4, 


+ .... + P p J ( l-P e )(!~) J + 


(3. 2. 2.1) 


or 


B x (s) = 




1 - pJ 


Pi 


e x s -h (i 1 


or 


B 1 (s) 


P 2 ( 1-P e ) 

s + 'p 1 ^ i- p~T 


( 3. 2. 2. 2) 


Thus, it is as if the service rate is deteriorated to, 


Pi eff = ^(1-Pg), the distribution 
exponential. V<e can now apply tho equations 
M/M/ l/N system. 


remaining 
for an 


3.2,3 Retransmission after Finite Time-out: 

Let the time-our = i . -~~™ as defined (Eqn. 3.1. 3) . 

Pi 

With time out modelled as Erlang-i process. Each retrans- 
mission, one service time (exponential) m series with an 
Erlang-i (i exponential servers m series) timeout stage , 
is effectively on Erlang~(i+l) stage. The L.T, of the 
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effective service time becomes (ref. Fig. 3.10 and eqn.3.14) 


B(s) = (l-p e ) ( , 


I»1 P, 

“ “> + Pe^-PeH^) 


i +2 


+ 


j, w ^ , j( i+l)+l 
+ ... + p e ( 1 -p) ( 7 ^ 7 ) + 


or 

B(s) 


(l-p e ) (ji x / s+p x ) 
[1-P e ( H 1 /s+H 1 ) 1+1 ] 


(3.2. 3.1) 


The effective service rate, is thus (ref. Appendix C), 


^leff 


1 

bTo} 




l+( i+l)(p e /l-p e ) 


(3. 2. 3. 2) 


We assume that the distribution may be assumed exponential. 
This approximation is valid because of the low value of p e 
under normal operation (Appendix C, eqn. C.18). 'We can then 
use as the service rate m the M/M/l/N model to evaluate 

the performance. 


3 . 3 SEND-AND-WAIT SCHEMES 


In these schemes the acknowledgement for a packet must be 
received prior to transmission of the next packet. This 
effect is incorporated by assuming no waiting room m front of 
stage 2 m Fig. 3.1. The effective service consists of stage 1 
server in series with stage 2 server ( pg = l)» with no waiting 
room in front of stage 2 server. 
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Let B^(s) represent the Laplace transform of the pdf 
of the effective service time at stage 1. This corresponds 
to either of the cases discussed ib Sec. 3.2, The L.T. of 
the pdf of overall service time, B(s) is given by, 

B(s) = B^s) (3.3.1) 

where the second term on the right-hand side represents the 
L.T. of the pdf of service time at stage 2 (exponential). 

The physical (successful) transmission from the node 
again corresponds to departure from stage 1. The service time 
at stage 2 is interpreted as delay for the next packet waiting 
for transmission. 

When retransmission are not present (p e = 0), the 
coefficient of variation of the effective service time at 
stage 1 is minimum (=1). The two exponential servers (stage 1 
and stage 2) in series gives rise to the coefficient of varia- 
tion of the overall service time less than 1. Therefore, 
modelling the system as M/M/l/N results in overprediction 
of the delay, nodal blocking probabilities etc. When 
retransmissions are present the hyperexponential nature of 
this stage brings the coefficient of variation of the overall 
service time closer to unity. M/M/l/N system analysis would 
therefore give better results. 
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3.4 CONTINUOUS TRANSMISSION WITH FINITE WINDOW AND NO 
ERROR RECOVERY 


Referring to Fig. 3.1 this corresponds to the case when 
P e « o. 

The effect of server 2 and the finite queues in front 
of server 1 and 2 prevents a closed form solution as in 
previous cases. A Numerical method is presented. 


We start with an initial value of pg e.g. P 2 = o. We 
incorporate the effect of blocking by server 2 by calculating 
the effective service rate of server at stage 1, as 


^leff 


H 1 P 2 


(refer Appendix C,eqn.C.10) 


( 3.4.1) 


For small P 2 values (less than 0.2) exponential distri- 
bution for the effective service time at stage 1 may be 
assumed (ref Appendix C, eqn. C. 14) . The traffic intensity 
at stage 1 is calculated as p ^ f Using the M/M/l/N 

equations we find p Q , the probability of server at stage 1 idle 
(no packets at stage l). The throughput from stage 1, X^, 
is given by, 


\ — P-^ e f|-( 1 -*Pq ) (3*4.2) 

The second stage blocks output from first stage with prob. 

P 2 « For X to be the system throughput with blocking 
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probability, the offered rate to the second stage, 

X 2 > must be such that, X = X 2 ( l-p 2 ) 

or 

x 2 = X 1 /(1-P 2 ) (3.4.3) 

Thus the traffic intensity at second stage P ^ may be 
calculated, 

P 2 = X 2 / ^2 (3.4.4) 

The new value of P 2 is calculated as the probability of W 
packets at stage 2, 


or. 


? 2 ( 1 - p ) 


P 2 


w+1 


(3.4.5) 


The equations (3.4.1) - (3.4.5) are iterated till convergence 
(which is possible if P 2 converges). The blocking probability 
at node 1, with the simplified assumption of Sec. 3.1, is 
the probability of N packets at stage 1, 


or. 



<=J( i- p ) 

1- PjN+1 



(3.4.6) 


Thus we have % 
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Algorithm 3.1 
begin 

Initialize = 0 
Until satisfied do 
begin 

p leff ^ 2 +p 2^1 
P 1 ~ X//p leff 

P 0 = ( l- p 1 )/C l— p 1 N+1 ) 

\ “ H lef f(l-P 0 ) 

2 = * lA i- P2^ 

P 2 ~ ^2/ ^2 

p 2 = P”(1-P 2 )/(1-P 2 W+1 ) 
end? 

P 1 = P l\ 

end? 

3.5 CONTINUOUS TRANSMISSION WITH FINITE WINDOW AND 

RETRANSMISSION WITH FINITE TIME OUT 

This is the most general case. We extend the algorithm 
developed m the previous section to account for the effective 
service at stage 1 with retransmission. This was dealt with 
in isolation, in Sec. 3.2.3. The effective service at stage 1 
is hyperexponential . tags* 2 


provides a server in 
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series, with a probability, p 2 * This has the effect of 
reducing the variance. Thus we expect the overall service 
time to be closer to exponential than the cases in Secs. 
3.2.3, 3.3 and 3.4. 


Vile also introduce, at this stage, the probability of 
redundant retransmissions, p . A redundant retransmission 
occurs where the acknowledgement is delayed beyond time out. 
This has to be accounted for in some of the previous cases 
we had deliberately omitted to avoid confusion. Appendix D 
gives the method for calculating p r (eqn. D.6). It is seen 
that p r is a function of p 2 which depends on p 2 » Thus 

algorithm 3.1 is to be modified to iterate over 
P r and p 2 * We also append the calculation of the parameters 
of interest to complete the algorithm. 


The effective service rate of server at stage 1 (not 
including blocking) is given using eqn. (3*2. 3. 2), 


^lloss 


l+( i+l)( p/l-p) 


where p = p + p . 


(3.5.1) 


p is constant but p changes m each iteration. Hence, 

^lloss Ganges in each iteration. The overall service rate at 
stage 1 including the blocking, Is found from servers in series 
consideration (eqn. 3.4.1), 



61 


, _ ^lloss ^2 

^leff “ — (3.5.2) 

^^lloss 

consideration of low p g and p 2 values lead to exponential dis- 
tribution approximation for the effective service time 
( eqn. C.16). Thus the general algorithm is ; 


Algorithm 3.2 


begin 

Initialize p r , p 2 = Oj 
While not satisfied do 
begin 


P = P e + P r 


P'lloss 


P 1 

i+( i+1) ( p/l-p) 


_ p lloss 1*2 
ff P 2 +p 2 |i Uoss 

h ~ ^ p leff 

P 0 = (i-pp/d-P^+i) 

X 1 = p leff (l -Po } 

X 2 = V<l-P 2 > 

Pg ^ ^ 2 / ^2 
p 2 = (i-p) p 2 w /( i-p 2 w+1 ) 
p r 


is given by eqn. (D.6) 
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end} 

P 1 * P l N p o 
Y - \ 

Pl [l-(M)P^ N P X N+1 ] 

Q = -1 l =. 

(i- pjMi- p^ 1 ) 


D = Q ± /y 

P = y/D 
end} 

where Y = throughput from physical node, 
- queue length at stage 1 


(eqn. 2.2.4) 

(eqn. 2.2.5) 


D = delay faced by packet at physical node. 
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CHAPTER 4 

DESIGN OF SIMULATION EXPERIMENT 

4.1 RUNLENGTH CALCULATIONS 

To validate the analytic model of Chapter 3, Monte Carlo 
simulation was used. The runlength was determined for 
95 percent confidence at 5.0 percent tolerance of the 
mean blocked waiting times. A relative tolerance value, 
rather than a fixed value, was used because of the wide varia- 
tions in waiting time with different arrival and service 
rates. 

The observations were divided into blocks to reduce 
the autocorrelation. The mean of the waiting time for block 
was taken as the observation for runlength calculation purposes. 
The number of blocks was determined using the criteria 
derived from the Central Limit Theorem [DEO 79.Chapt. 7], 

It was observed that block sizes of the order of 200 were 
necessary to ensure little correlation between the mean of 
waiting time per block. 

4.2 IMPLEMENTATION OF SIMULATION MODEL IN PASCAL 

The names of variable used in this chapter are based on 
the notations of Chapter 3. An array INPUTQUEUE , of size N 
(physical buffer size) stores information on packets waiting, 
being transmitted (in service) or waiting for acknowledgement. 
The information stored is the time at which a packet is to be 
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packet is to be appended to the table (This will be explained 
in the next few paragraphs)* 

A routine INSERTQUEUE accepts the array on which 
updating is to be done, the time of arrival the maximum 
number of entries as parameters* It updates the table by 
finding the position j, such that ine last j packets have 
departed prior to the arrival. It calls a routine 
REMOVEQUE. This routine accepts as parameters the array 
on which updating is to be done and the number to be 
removed in a FIFO manner. This is done by moving the entries 
j places up starting from the (j+l)st position. It 
returns the number currently in the queue. This is sub- 
sequently returned by INSERTQUEUE a calling sequence. 
INSERTQUEUE also returns the position where the information 
on the current arrival is to be placed. INSERTQUEUE signals 
via a boolean variable QSTATUS whether the packet can be 
accepted or not, due to buffer limitation. 

The waiting time for transmission is calculated 
recursively. Assuming the current packet under observation is 
(n+l)th packet, the waiting time, 

W n+1 = max(0, TC n - a n+1 ) (4.2.1) 

where TC n corresponds to time of successful transmission of 
the previous packet and a n+1 the time of arrival of current 
packet. 
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The (n+l)th packet is ready for transmission at 

time, 


y n+l a n+l + W n+1 


(4.2.2) 


The 2nd table ( array) , ACKqueue, of size W stores the 
time of receipt of acknowledgement, (departure time from 
node) for packets waiting for acknowledgement at last 
update. This is scanned to check if transmission of the 
packet ready for transmission is to be blocked because of 
window constraints. The same two routines INSERTQUEUE and 
REMOVEQUEUE are used. The parameters passed now are the 
array ACKqueue, the time of arrival as Y n+ ^ and the maxi- 
mum length as W. The queue length ( Q 2 ) is returned. 


If W packets are occupying ACKqueue, the time of trans- 
mission attempt of the packet ready for transmission is 
delayed by a factor <5 n+ -j_ accounting for the blocking. This 
is accounted m the waiting time, by incrementing it 

by w 

Or, 


W 


n+1 


W n+1 + 6 n+1 


(4.2.3) 


where 6 n+ ^ = time of receipt of acknowledgement of the first 
packet in ACKqueue - Y n+ -j_ (4.2.4) 


The first packet in ACKqueue is removed by calling procedure 
REMOVEQUEUE. 

The time of start of transmission for the packet ready 
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is modified to, 

Y n +1 = *n+l + 6 n+l (4 - 2 * 5) 

The transmission process starts now. The time of 
transmission (service), S n+ ^, is determined by generating 
a random variable from exponential distribution with 
parameter p,^. It may be noted that = pC where l/p, is 
the mean packet length in bits and C, the transmission 
capacity of line bits per second. 

The transmission process checks if a retransmission 
due to physical errors m transmission of packet, the 
acknowledgement or blocking due buffers being full at 
next node is necessary. Further more after a successful 
transmission and acceptance, several redundant retrans- 
missions may take place because of the delay in receipt 
of the acknowledgement beyond time out period. 

A procedure, RETRANSMIT, accounts for all of these 
effects. The time out factor, i, specified by the user, 
is used to calculate the constant time-out period. 

Time-out = — • i = — (4.2.6) 

The transmission of a packet, attempting transmission 
at Y n+1 > is completed at time, 

TC n+l = Y n+1 + S n+1 


(4.2.7) 
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A random variable is sampled from uniform distri- 
bution in (0,1) to determine whether the transmission, 
acknowledgement contained errors or the packet was not 
accepted. This is determined by checking whether the 
sampled value is less than the probability of the above 
events, specified by user. 

If the transmission is deemed unsuccessful, the 
next transmission is attempted at time, 

Y n+1 = TC n+1 + time out (4.2.8) 

This process is continued till successful transmission. 

The calculations on the time of receipt of acknowledgement 
is now done. This also includes the effects of redundant 
transmission. 

The acknowledgement for the packets are assumed to be 
received in the same sequence as their transmission. The 
acknowledgement intervals are sampled from an exponential 
distribution with the parameter, 

The time of receipt of acknowledgement, 

4+i ■ ma *< D n 2)+S nri' TC n + l + < 4 - 2 - 9 > 


where 
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( 2 ) 

' = time of receipt of acknowledgement for the 

previous (nth) packet, 

S>n +1 = the acknowledgement interval for the ( n+l)th 

packet. 

( 2 ) 

If is greater than the next retry time for 

the packet, y n+ ^, another retransmission is performed. 
Transmission is complete at, 

TC n+l = Y n+1 + S n+1 
The next retry time, 

Y n+1 = T( ^n+1 + ' t ^- me ou ' t calculated. 

( 2 ) 

If this is still less than , another retransmission 

is attempted. This process continues till the next retry 

( 2 ) 

time Is greater than These retransmission represent 

f 2) 

the redundant transmissions. is the entry stored in 

the tables, Inputqueue and ACKqueue for the current 
packet. All the acknowledgements received during redundant 
retransmissions are ignored for the purposes of the 
packet(n+lth) under observation. 

For the next packet however, the recursive equation 
(4.2,9) assumes the time of receipt of acknowledgement of 
the previous packet as the last of the acknowledgements 
received. We could have made more realistic assumptions. 
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However, the effect of these retransmissions means 
delay for subsequent packets. The effect of the window, 
postpones the delay to the (n+l+W)th packet. Since the 
statistical behaviour is of interest, incorporating the 
delay to the next packet, makes no difference. 

4.3 RANDOM NUMBER SEQUENCES USED 

Five different uniform random number sequences, 
borrowed from FORTRAN library, were used : 

a) RNDY1 - this is initially seeded by the user 
interactively. Samples from these random variables are 
used to seed the other sequences. 

b) RNDY2 - to generate interarrival times of packets 

from exponential distribution with parameter . 

c) RNDY3 - to generate transmission time (service time) 

of packets from exponential distribution with parameter 

d) FNDY4 - to check for loss of packet in transmission, 

loss of acknowledgement or blocking at next node,. 

e) RNDY5 - to generate the time to transmit the acknow- 

ledgement from adjacent node 

This is sampled from an exponential distribution with 
parameter j^* 
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4.4 DIFFERENCES BETWEEN SIMULATION AND ANALYTIC MODEL 

Simulation models must be kept as close to reality as 
possible. Assumptions in analytical models, however, have 
to consider the ease of mathematical tractibility. 

In the present study the points of similarity between 
the analytical and simulation models are i 

i) the mterarrival times of packets are independently 
distributed random variables from exponential distri- 
bution, 

ii) the transmission times (based on packet lengths) are 
independently distributed random variables from 
exponential distribution. 

iii) the time for acknowledgements are independently 
distributed random variables. This is true ev^en for 
retransmission. The last assumption in the simulation 

model is based on the observation that the acknowledge- 
ment is strictly not governed by the length of acknow- 
ledgement packet. Other factors like delay in proce- 
ssing at next node, transmission and propagation delays 
and piggybacking of acknowledgements are also involved. 

The points of difference are : 

i) The simulation model uses a constant time out period 

while the analytical model assumes the time out to be a 
random variable sampled from a Erlang-i distribution, 
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where i is the time out factor. 

ii) During a retransmission, the originally sampled service 
time (time for transmission of packet) is used in the 
simulation model. This is because the physical packet 
length cannot change. The analytical model assumes 

an independently sampled value of service time 
(stage l) for each retransmission. 

iii) Dynamic sharing of buffers between stage 1 and stage 2 
is assumed in the simulation model. The analytical 
model assumes completely partitioned buffers where N is 
available from incoming traffic and W for packets 
waiting for acknowledgement. This is based on physical 
buffers size, N being sufficiently large as compared to 
the expected occupancy by unacknowledged packets, in 
network nodes. 

iv) In the simulation model, the actual blocking time is 
calculated for a packet that cannot be transmitted 
till the receipt of the acknowledgement of the first 
packet in the acknowledgement queue, when the number 

of unacknowledged packets is equal to the window size. 
In the analytical model, the Markovian assumption on 
the acknowledgement intervals is made. This means the 
blocking experienced by a blocked packet is on the 
average, mean acknowledgement interval. This is because 
the residual service time of the first packet at stage 2 
is independent of the time already spent in service. 
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CHAPTER 5 

DISCUSSION OF RESULTS 

5.1 ACCURACY OF THE MODEL 

The results obtained by simulation and analysis 
(Alg. 3.2) are listed in Table 5.1* Results of a second 
method of analysis, which incorporates an approximation to 
account for the hyperexponential nature of the server at 
stage 1 are also listed. Fig. 5.1 clearly indicates the 
need for the second method under moderate traffic. This is 
further supported by results of Appendix B (eqn. B.6). The 
approximation is explained in Sec. 5.2. 

The analytical methods are able to predict reasonably 
well the parameters of interest at a computer network node, 
namely, throughput, delay, network power. As explained the 
second method gives better results. Also the buffer overflow 
probabilities are approximated reasonably well by the second 
method. The first method is not so accurate for this para- 
meter. It may be seen from Figures 5.1 - 5.5 that in most 
cases the two analytical methods outline the envelope m 
which the actual values, given by simulation, lie. 

The mean and standard deviation of the effective service 
times ( x , a ) and interdeparture times (d, cr d ) at stage 1 
are predicted fairly accurately, in most cases. * <=Vu e ff) 
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and d (= l/y) are given by algorithm 3.2, a is found 
using the coefficient of variation, C^, of the service time, 
which is given by eqn. C.18 (Appendix C) .'Ve have omitted the 
effect of P 2 (eqn. C.18). This becomes appreciable only 
for W in the range 2-3. 

The coefficient of variation of the interdeparture 
times, C d » is used to determine c^. is approximated by 

the result for the interdeparture times m an M/G/i/“ system 
[MARS 68], 

2 ^ 

C d = 1+ P^(C b -l) (5.1.1) 

where 

P = X x 

5.2 APPROXIMATIONS TO ACCOUNT FOR HY PE REXPONEN T I .\L SERVER 

The stage 1 server, effectively being hyperexponential 
(Appendix C, eqn. C.16) makes the 2 stage model 
M/H 2 /l/N + . /M/l/W, with successful transmission of 
packets, corresponding to departures from stage 1. Two 
approximations are used s 

i) The second stage is effectively an Hg/M/l/W system, 
due to output of stage 1 being hyperexponential 
(eqn. 5.1.1). The coefficient of variation of the 
arrival process to second stage is used to calculate 
the parameters of the H 2 distribution (eqns. A.l, A. 5). 
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The traffic intensity, cr, is found by eqn. B.4. The 
M/M/ l/'U system analysis is used to calculate the queue 
length with this traffic intensity, a. The approximation 
made is the increase m queue length due to greater 
variation m service time is incorporated as increase 
in traffic intensity m an M/m/i/W sys tern. 

li) The hyper exponential server at stage 1 gives rise to 
waiting times greater than that for an exponential 
server with the same mean* The first stage is effectively 
M/H 2 / 1 /N. However, H/G/l analysis is more difficult than 
Gl/M/l analysis as pointed out in Sac. 2.3. V ! e incorporate 
the effect of increase m queue length duo to the 
increased variation in service tome as if the M/M/l/N 
system was having a higher traffic intensity. >/i/e use the 
h 2 /m/i/« correction again (eqn. D.16) again. To do this, 
we reflect the coefficient of variation of the service time 
at stage 1 to the arrival process and assume the service 
time is exponentially distributed. 

Considerable approximations have been made. These require 
careful scrutiny. A search in literature, or analysis, for 
queue lengths for M/H 2 / 1 /K and Hp/M/ 1 /k systems, and coeffici-- 
ent of variation of interdeparture times from anM/H 2 /l/K 
may be done for better results. 
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5.3 SENSITIVITY OF PO-ER TO VARIATIONS IN PREDICTED 
VALUES OF PARAMETERS 


is, 


The network power, referring to notations of Chapter 3, 


Throughput 
' delay 


or 



( 5.2.1) 


P is therefore more sensitive to errors m the prediction 
of y. It may seen that y = l/’d, is predicted reasonably 
accurately in most cases. The error in prediction of P, in 
most cases, is due to the incorrect prediction of L. 

i) This may due to the assumption of buffer occupancy by 
unacknowledged packets small enough to be neglected 
(Sec. 3.1). The error introduced by this is not expected 
to be substantial, simulation results show the expected 
occupancy by unacknowledged packets Q 2 ? do not exceed 3. 

With buffer sizes of the order of 30 or more the error 
in N for queue length calculation (eqn. 2.2.4) is not more than 
10 percent. Further more the analytical methods underestimate 
Q 2 « This appears to be due the approximations used. 

ii) The mam cause for discrepancies in calculation of is 
due to the approximations made to incorporate the effect 
of H 2 server at stage 1. 
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It may be noted that the correction factor for the 
hyperexponential server is not incorporated m the second method 
for traffic mtensties at stage 1, > 1. This is because the 

result derived is for H^/M/l/® system, where the ergodicity 
condition is < 1. Studies may be done to see whether 
the correction factor is significant for ^ > 1. We do 
not expect this to be significant due to heavy traffic. At 
= 1, the correction factor is zero (Appendix B). 

5.4 DESIGN CHOICE OF BUFFER SIZE AND WINDOW 

The second method is found to be more accurate in 
predicting the buffer overflow probability, p^. This may be 
used to find a proper buffer size, N. 

A plot of p^ y s N is shown in figure 5.2. Plot of power 
(P) vs N would give false notion of increase m power with 
reduction in N, This is due to reduced queue lengths and 
hence delays. A node however cannot be considered in 
isolation. Smaller N values give rise to higher p^ values. 

This means increase in retransmissions and hence delay, at 
previous node. 

The choice of N should be one that reduces buffer over- 
flow probability within acceptable limits. At the same time,N sho- 
uld notbe too large. The response to congestion would be 
sluggish. It may be noted (Chapter l) that information (effect) 
of congestion in an interior node must be passed to the 
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entry points as quickly as possible for good network flow 
control scheme. 

Variation of power vs W is shown in Fig. 5.3. There is 
no optimal W. However, beyond a certain range, it is seen 
that the increase in power with W is very small. Vv must be 
chosen just sufficiently large to give the desired performance. 

A substantially large W increases the probability of congestion. 

Furthermore the model allows a calculation of the upper 
bound of the blocking probability, p^ due to window. This 
happens when the traffic intensity (note traffic is already 
regulated by stage 1 and normally stage 2 server is faster 
than stage 1 server) at stage 2, P 2 is 1. p B is then . 

For W = 8 (.e.g. HDLC), p B = l/9. This is strictly an upper 
bound. P 2 is actually much less than 1. 

It is further observed that the values of N and W on the 
performance are correlated. A small window would restrict the 
throughput irrespective of the buffer size. It is observed 
from the results that beyond a certain limit e.g. C'i = 5 
increase of window increases power as 1 a fine tuning mechanism' . 
The choice of W = 8 for HDLC is thus a good one. 

The design procedure for N and W is thus : 

Given traffic parameters, time out, retransmission 


probabilities . 
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i) Initialize W between 5 to 8. 

ii) Increase N till approximately the acceptable 
buffer over-flow probability is attained. 

iii) Increase W for better performance, at the same time 
keep within reasonable limits for congestion 
avoidance. 

5.5 OPTIMAL TIME OUT FACTOR 

The analytical model predicts the performance with 
various timeout intervals close to simulation. Some cases 
are plotted in Fig. 5.4. A small time out will give rise to 
greater probability of retransmission, A large time out 

will degrade the performance of the transmission stage. 

Referring to sec, 3.5., 

111 eff i+(i+D 

where 

P = P@ + P r 

given by Appendix D (eqn. D.6), is a strictly 
decreasing function of i. Th: optimal value of i such that 
H eff is maximized may bo founo analytically or graphically. 

It was observed that for W > 10, approximation for 
W = 00 (Appendix D) may be used. 
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5.6 INDEPENDENCE OF OUTPUTS 

It was observed that the autocorrelation coefficients 
( lag > 1) for departures from the transmission stage 

(stage 1) were sufficiently small. This allows assumption 
of independence of the outputs from engineering view- 
points. This is useful for extension of analysis to network 
of queues. The inputs to the next node may thus be assumed 
independent. 

This is an important observation because of Finch's 
result [FINC 59] that the output of an M/G/l/K system is 
independent only if the server is exponential and the waiting 
room, K = <* , We may look into the method of proof for 
correlation of interdeparture times in an M/M/l/N system 
to get an insight to the observed results. 

Finch showed that, for the output to be independent, 
the probability of a departure leaving no customers behind, 

p 0 = 1- X X ( 5.6.1) 

’h , X being the notations introduced earlier. However, 
M/M/l/N system analysis leads to 

P 0 = (1- P)/(l- P N+1 ) 

= Xx . 


( 5 . 6 . 2 ) 


where P 
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Thxs violates eqn. (5.6,1) and hence the original assum- 
ption of independence was false. However, from engineering 
considerations, for large N ( > 20) the denominator in eqn, 

( 5,6,2) is sufficiently close to unity, particularly if f is 
not close to unity. 

Furthermore, the deviation of coefficient of variation 
of the server at stage 1 is not substantially different from 
unity, as shown by the results. 

Also the coefficient of variation of the interdeparture 
times is bounded by the coefficient of variation of the 
service time, if the arrivals are poisson (eqn. 5.1.1). The 
upper bound occurs only when the server is busy at all times 
(X = “). Thus independent arrivals with poisson distribu- 
tion may be assumed for the next node. 


Further work on the autocorrelation coefficients may 
done based e.g. on [DALE 68] or Cox’s result [BURK 72] 


1 

y=l 


p 2 ( i-c 2 ) 

2[l-P 2 ( 1~C 2 ) 


( 5.6.3) 


where a r = autocorrelation coefficient of lag r 

C = coefficient of variation of service time 
P = X x 
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5.7 REPLACEMENT OF EFFECTIVE TRANSMISSION PROCESS BY A 
HYPEREXPONENTIAL SERVER 

The transmission stage (server 1} with retransmissions 
after time out is modelled by an equivalent hyperexponential 
server (eqns. C.16, A. 4). 

Simulations were carried out, replacing the transmission 
stage including retransmissions after time out mechanism, by 
the equivalent hyperexponential server. The mean and variance 
of service time of the server was fed from the results of 
simulations with the retransmission mechanism. It was observed 
that the results were accurate to those of the original model. 
The blocking probability by server 2 is slightly higher with 
the replacement by hyperexponential server. 

This has two consequences : 

i) Any complex transmission/retransmission mechanism may be 
given an equivalent representation. Queuing systems 
analysis require only the mean and variance. Approximating 
to the n moments requires an H n representation. Hence, 
the H 2 representation is sufficient. 

ii ) Considerable saving m computer time results by replacing 
the complex retransmission procedures m simulations, 
particularly when simulating network of queues. 
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This however requires a good approximation for the mean 
and variance of each transmission stage in the network. This 
may be taken up for further study. 

Simulations wore also performed under Erlang-3 assumptions 
on packet transmission times. This is based on physical 
packet length having less variation than that given by the 
exponential distribution. This stage, with coefficient of 
variation < 1, in series with the retransmission stage, 
coefficient of variation > 1, gives overall performance closer 

to exponential. The M/M/l/N** */M/l/W analysis was found 

more accurate in this case. 

Simulations were also performed varying x , (i^ and p 2 
but keeping the ratios fixed. The parameters which lead to 
calculations of delay, throughput etc., namely, pl,p2,Ql,Q2,pr 
were quite close. This was performed to check for the effect 
of eff - (sec 3.5). This gives rise to P 1 = X / not 

strictly dependent on the the ratios. The explanation is 
that for most cases p 2 is sufficiently small. It is expected 
that for p 2 large (e.g. W = 2-3) change of the parameters x 1? 
p 1 ,}i 2 with ratios constant may give results that are not very 
close. 

Under operating conditions of W > 5, the results m 
Table 5.1 may therefore, be used to infer performance measures 
for other ranges of values of X 
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5.8 STABILITY OF THE NUMERICAL METHOD 


It was observed that, for the cases considered, conver- 
gence was independent of the starting values of p 2 and p . It 
was found that m some cases, however, convergence was not 
possible because of oscillations m the values of p . It was 
observed, in most of these cases that the mean of the two 
values between which p r was oscillating, was close to the 
actual value, given by simulation. These cases are listed 
below : 


l) * = 7.0, p 1 = 10.0, p 2 = 10.0, N = 15, W =10, i = 5 
for values of p e < 0.05. 


ii) * = 40.0, = 50.0, p 2 = 50.0, p 0 

i = 5. 


0.01, N = 15, W = 10, 


5.9 SOME OTHER POSSIBLE AREAS FOR FURTHER WORK 

An ar.-a which needs attention is the modelling of the 
acknowledge process (stage 2). With piggy backed ACK schemes, 
and acknowledging upto K previous packets 

e.g. HDLC, it is felt that exponential server with bulk service 
[KLEI 75] would be a better model. 

Analysis for state dependent processes can be done by incor- 
poration m the current numerical procedure. These can be useful 
example determining loss probability, of voice packets where the 
probability of arrival of a packet and duration of service may be 
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state dependent on the number of packets already present 
(no. of active talkers). 

Extension of dynamic flow control measures based on 
Markovian Decision model of Kemani and Kleinrock [KERM 80a] 
may perhaps be done upto link level. It may be interesting 
to see whether a dynamic control on link level window, time outs 
would actually improve the performance as in end-to-end case. 

The effective transmission time may be interpreted due to 
increased average packet lengths and no retransmission. A 
comparison with error coding methods e.g. forward-error- 
correction for random or burst errors, may be done for noisy 
channels . 

Finally, the extension to the overall performance of a 
network needs to be done. Direct application of Jackson’s 
theorem, BCMP queues are not applicable due to finite queuing 
assumptions. Without this assumption the effects of retrans- 
missions are lost. 

A generalized numerical procedure extending the present 
one ofor a single link may be designed. 

For example for a tandem network we would have to modify 
algorithm 3.2 to an algorithm that would be an extension of 
the method of Irland et al [IRLA 77], 



86 


Now there are two levels of convergence. The outer 
loop would account for the nodal blocking probabilities 
for the tandem physical nodes. For n nodes start with 
Pl = 0, j = 1, . .., n. Pij corresponds to Pi determined 
by algorithm 3.2* for node j. 


If the arrival rate to the network , then the 

effective arrival rates to individual nodes get modified 
due to retransmission to [IRLA 77] 


X l“ 


X = 

3 


( 1 - 


* 11-1 
■PiP 


2, . . . , n 


Using this value of for ^ m Algorithm 3.2, we now 

iterate for stability for each node. Note that at the 
beginning of the iteration of the inner loop (Algorithm 3.2) 
for individual node. 


P 


ej 


= Ft. + P 


1J 


where 

rate. 


p T is the physical packet or acknowledgement error 
L J 

on link out of node j. 


Thus, we propose the following numerical procedure for 
tandem network. 


Algorithm 5.3 

Given a network of n nodes in tandem with arrival rate 
X, time out factors i., 

line loss rates p^ , transmission rates acknowledgement 
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rates buffer sizes n , window 

^ j * 3 ” 1,2, « • • , n . 


Initialize = 0, j = i, . n. 

Repeat 

For j=l,2, .. ,n do begin 
If j=l \ olse V, = Xi^.i/U-Pu 5 

(* Algorithm 3.2 *) 


Initialize p^ _ = p 2 ^ =0, j = 1, n. 


P ej = 
While 
begin 


P l3 + P 13* 
not stable p 

r 3 



do 


P J = 


ej 


+ P 


ri 


^HOSS,3 


l+( ij+ 1 ) 



p lef f , j 


p lloss.i ^2 j 
^ 23 +p 23 ^lloss ,j 


P l3 P 3 //fp lef fj 

N +1 

P oj = (1- P l3 )/(1- ) 


*lj ^leff , 3 ^"~ P o 3 ^ 

X 23 = X 13 /( 1 “ P 23 } * P 2j = X 2j /p 2j 


sizes 
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end} 


w w +1 

p 2J = d- p 2j ) p 2 //u-p 2 ] ) 

Calculate p r ^ by eqn. D.6 (Appendix D). 
end 

. N 

P lJ = P o 3 B 1 

until Pj_j* s stabilize 


for 3 - 1, . . . , n 


begin 

calculate by eqn. (2.2.4) using P^,N^ 
Dj = where M = X ^ s: system 


end 



89 


5.10 SUMMARY 

Performance evaluation of existing flow control schemes 
or choice of parameters for optimal performance for a new 
scheme are required in computer networks. Simulation offers 
one solution but requires large amount of computer time. Exact 
methods, based on queuing theory are complex for computation 
purposes. Furthermore, the effort spent m exact solutions 
may be futile from engineering viewpoints due to the discrete 
values of physical components e.g. transmission line capa- 
cities, timeout intervals etc. The need arises, therefore, 
for simpler approximate methods • 

The approximation method (based on a 2-stage queuing model), 
presented in this thesis, is able to predict the performance 
parameters reasonable well. Furthermore, each control 
parameter is given an interpretation in terms of queuing 
effects e.g. retransmission after finite time out increases 
the coefficient of variation of the transmission stage, the 
acknowledgement process reduces the coefficient of variation 
of the transmission stage, window limits the waiting room m 
front of the second stage etc. 
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APPENDIX A 

DETERMINATION OF PARAMETERS FOR AN H 2 
DISTRIBUTION GIVEN MEAN AND COEFFICIENT 
OF VARIATION 

Let the mean be* x = ' l/\x 

and coeff. of variation = C^ 

The H 2 system, m Fig. A»l, is represented by two exponent 
tial processes m parallel, with a, (l-a) the probability 
of service at the respective server. 

t 


x 

Fig. A. 1 

Then [WHIT 75, pp 153], 

= 2a p, P 2 = 2 (l“ a ) P (A.l) 

Now, the expectation, 

Ei = ct • 1 + ( 1—cc ) _1 

P 2 



1 

P 


(A. 2) 
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The 2nd moment is given by [KLEI 75], 


d 2 \ + ( 1-a ) \x r 

ds i orri) Ti hCT 


L 


2a 

M-i 2 


Zllzzl 

\x 2 2 


s=0 


( A. 3) 


The variance, 
V. 


2a + 2 ( 1-a ) - E 2 






= 2a 


2 2 


+ 2 ( 1-a ) - 1 


4 ( 1-a Vh? 




=1 [ ( 1-a ) + a 


|A 


2 2a (1-a )' 2a Tl-a) 


1 1 


= 1 [ 1 - 1 ] 

..2 2a"~CT-a7 

{J, 

or , V t = E 2 [ l-2a ( 1-a). ] 

z Z Sad^aT^ 


or, C 2 = V. = lzMl=Zl 

D — o 2aTI-a) 


(A. 4) 


To solve a m terms of we have 


2a ( 1-a) C b 2 + 2a ( 1-a) -1 = 0 



or, 2a 2 [ 1 + C b 2 ] ~ 2a [l+ C b 2 ] +1=0 
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or, a - a + 1 1 


2 rr+ c 2 3 * o 


This is solved to give, 

a = i c 1 ± h 


= i [ i ± /c k 2 -i 


2(1+ C b ) 

] 


v C b 2 + 1 

Choosing the -ve sign to keep a <0,5 [WHIT 75] 
we have, 


a = 1 r i 

■55 1- 1 


c b 2 - 1 
C b 2 + 1 


] 


(A. 5) 


Queuing system analysis requires only the mean and 

variance. Thus any system with coeff, of variation, 
c b > 1 can be given an H 2 representation by eqn. A. 5. 
If approximation to n moments were necessary, an equi*- 
valent H n system would have been required. 
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APPENDIX B 

TRAFFIC IN1ENSITY FOR MI H^-M/ j_/«> SYSTEM 


Let the 

arrival 

process, assumed H 2 , have the 


transform. 

A*(s) 

= “JH _ + '(l-a) p 2 

(s + \x ± ) Is + p 2 f 

(B.l) 

where 

^1 

= 2a A =2akp 



P 2 

= 2( l-a ) *= 2 ( l-a) kp 

(B.2) 


(Refer to Appendix A). 

where l/ \ is the mean mterarnval time, the service 
time follows exponential distribution with mean l/p. and 
k = ^/m-* k corresponds to the traffic intensity of M/M/l 
system if arrival was poisson. 


A* (s) = s [ a|i 1 + (l-a) p 2 ] + p 2 

(s + (s + p 2 ) (B.3) 

Now for a Gl/M/l system the traffic intensity, a is 
such that (eq. 2.3.2. 4) 

a = A* ( \x - per) 


This gives the state probabilities, 
p k = 0 k (1 -o k ) 


, k as 0, 1, 2, 


• « • * 
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a - \i ( 1 - 0 ) [2a 2 k\x + 2 (l-a) 2 k|i] + 4 (l-a)ak 2 p 2 
[p (l-a) + 2a kji] { ji( l-a) + 2 ( l-a) k [a] 

2 

Cancelling pi and letting r = ( 1 - 0 ), 

r = 2 rk [a 2 + (l-a) 2 ] + 4 k 2 a (l-a) 

[r + 2ak] [r + 2 (l-a) k] 

= 2 rk [ a 2 +(l-a) 2 ] + 4 k 2 a (l-a) 

r 2 + 2 kr + 4a (l-a) k 2 


Rearranging and cancelling terms, 

r 2 - r 2 - 2kr 2 - 4 ra ( l-a ) k 2 + 4 rka ( l-a ) = 0 
Cancelling r throughout and rearranging, 
r 2 - [l -2k] r - 4a (l-a) k (l-k) = 0 

r = [l~2k] + i (l-2k) 2 +[l6 k (l-k) a (l-a) ] 

2 


reverting back to 0 , 

0 = 1 + k t 1 f ( 1-2 k) 2 + 16 k ( 1 - k),a ( l-a) 

2 2 

We choose tho - ve sign so that 0^1 

or, 0 = k + [ 1-1 / (1-2 k) 2 + 16 k (l-k) a (l-a) ] 

2 2 

(B.4) 


= k + H 2 c 

where H 2 c = correction term due to H 2 input, k = 
the traffic intensity for M/M/l queue. 



Thus given the coeff. of variation of a process, 
we calculate [Appendix A] 
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a 


1-1 
2 2 


' C ba 2 


C ba + 1 


and then H 2 c = I [l - /( 1 - 2 K) 2 + 16 k ( 1 


when a = 0.5 


k) a ( 1 - a)] 
(B.5) 


'/(l-2k) 2 + l6k ( 1 - k ) a ( 1 -a ) 
={(l - 2 k) 2 + 4 k (1 - k) 

= / l-4k + 4k 2 + 4k-4k 2 =l 


or, H 0 c = 1-1=0, (as expected - this is the M/M/l case). 

2 7k 2 
when a = 0 

H 0 c = 1 [l - V( 1 - 2 k) 2 ] = 1 [l - (1 - 2 k)] 
z 2 2 


= k 

Thus the upper bound for traffic intensity for a H 2 /M/l/ «, 
queue is twice that for M/M/l/<» queue with some parameters. 

Also, 

H 0 c = 1 [l - { (l - 2 k) 2 + 16a (l - a) k (l - k)] 

2 2 

for a constant a, let a (l - a) = b 

H 2 c = 1 - V 1 - 4 k (T^Tb) + Tk 7 (1 -4b )] 

= 1 [1 - f l - 4 k (1 - k) (1 > 4 bj J ( B .6) 

2 

This is maximum at k = 0.5. 

Thus the effect of this factor is more pronounced at 


moderate traffic. 
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MEAN AND VARIANCE OF EFFECTIVE SERVICE TIME 
C.l EQUIVALENT REPRESENTATION 


In this esection we derive the expression for the 
mean and variances of the effective service time at stage l 
of the model m chapter 3. 

Referring to Fig. 3.1, each retransmission means an exponen - 
tial stage with parameter p^ m series with a timeout stage 
which is Erlang=i (i exponential serves with parameter p^) 
in series. Effectively the retransmission s ta e is there- 
fore Erlang - (i + l). Letting j = i +1, and using nota- 
tions of chapter 3 (ref. se c 3.5 ) the overall 

system is shown m Fig c.l, 


j — (^£*'La.r\q -L )t 


— 
& V/i-SfG/v/ 


hp*. 


-> 


r 

ACK 


Fig. c.l 

p is the probability of retransmission 0 
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To allow calculations of the mean and variances of 
retransmissions only, overall service with no retrans- 
missions (See. 3.4) wj modify tho representation tos 



fig. C.2 

This docs no L change tho moments of the overall service 
time pdf. he shall uso tho fact that the L-T. of pdf of 
servers m series is the product of the L-T.'s of indi- 
vidual servers. This loads to the variance of two 
servers m series is the sum of the individual variances. 
'A'o shall calculate the overall variance by calculating 
variances of tho individual stages and B 2 . 

The L.T. of the retransmission stage, 

(s) = (C.l) 

1 - p (ll_ ) j 

S + 

Taking 1st and 2nd derivatives, 

(s) = - (i - p ) j P p x j ( s )j+1 

[ i - p ( — 3 2 

s + 


(C.2) 
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B r <«> 


\ - [1 ~ p ( — c( i~p)pj( j+i) 




*1 V S + \1. 


j + 2 
) ] 


(l-pjjpty 1 2[l-p ( __^i ^ ] x 


-P( 


S + Ph 




[i-P (sTi ^-) J ] 4 


( C. 3) 


The moan, 


«1 = “ B 1 (°) = llzRhs, =j£_ 


P 1 ( 1-p) (C.4) 


The 2nd moment, given by B ”( o ) } 

J ~s? = IP (l+.l+.ip-p) 

P-j 2 (1-p) 2 

The variance, 

V( s x ) = ^-^ 2 
(1-p) 


(C.5) 


(C.6) 
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The coefficient of variation, 


c b ( s i> =VWi7T~ = 


+( i-p) 
jp 


(C.7? 


The general nature of *' > C fa for the retransmission stage 

is shown below; 

/ 

C b 

1,0 

o p TTo 



The second stage gives rise to, 

P x ( 1-P 2 ) 




B 2 (s) 


( s +P-i ) 


p 2 ( 


•) ( 




p 2 


s+p_2 


•) ( C.8) 


W 


I t 


e calculate B 2 '(s) and B ? (s) by differentiation of 


o qn . C . 7 


(C.9) 


The mean, 


-B 0 (0) = _1_ + p 2 

Pi p 2 


(C.10) 


The 2nd moment 


s 2 = B 2 = 


2 . 


£t + 2p 2 T ^2 + 2p 2 <C ' U) 
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from oqn C.10, we have , 


1 

1 ->■* ^2" 


1 

^2 


2 + 2p 2 7TX, + 


the variance, from eqns. (C. 11) and (C.12) 


(C. 12) 


V(s 2 ) 




- s, 


~S + 2p 2 ~2 
p l p 2 


2 1 

P 2 - 2 
P 2 


"^r + P2(2-P2) ~2 

Pi P 2 


(C. 13) 


This gives, 

2 

c b ( s 2 ) = 


rj ^ 

~ — t~T 

P 2 +2 p 2 PiP 2 +P2Pi 


(C.14) 


When |i^ = p 2 = P» 


C b ( S 2> 


yi+p 2 ( 2~p 2 ) 

( h-p 2 ) ' 


(C.15) 


The cooff. of variation decreases from 1.0 to l/V2,as shown, 



p 2 1.0 


O’ 
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The parameters of fee overall service time at stage 1 of 
model are given by considering those of B and B, in 
series* 

The mean 

3*8, + S 2 
or 

®=ii^) + ^ + P2 ^ (c - ,6) 

The variance 


V(s) « V(s,) + V(s 2 ) 


_ 3P.( H.i-P) , 

~ O O * 

1^(1 -P) 



P 2 (2-P 2 ) 



Therefore, the coefficient of variation, C b (s) = 
The case p =0 corresponds to no retransmission. 


(G. IT) 

JWb) 

mm • 

s 


This is same as eqn. (C.14). 

% 

To find the coefficient of variation of server with 
retransmission (stage 1 of model in Chapter 3), we set 
P 2 = 0 in eqns. (C.l6), and (C.17). 


Therefore , 


& 


a+i-p)+0-p)‘ 


3P + ( i-p) 


(C. 18 ) 
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The general nature of the curves is plotted below. 



This shown that the server's coefficient of variation 
deviates from 1.0 at a certain range of p values. Thus 
physical network nodes with low retransmission probabilities 

may bo modelled as exponential. 
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APPENDIX D 

THE PROBABILITY OF REDUNDANT TRANSMISSIONS 

We base this calculation on the model and the 

notations of Chapter 3. 

The probability of redundant transmission, 

¥ 

P r " 2 p r ob (<K acknowledgements m Timeout) x 

k=1 

prob(K packets waiting for ack. ) 

(D.1) 


or, 

W-1 

P r = 2 prob(K ack's in Time-out) x prob 

k= ° ( >K waiting for ack.) (D.2) 

The last. form is preferred because of ease of 

calculation. 

Probability of K ack's m time-out, 

time-out 

p _ • - — ( \i 9 time-out) (D.3) 

because the acknowledgement process is poisson with 
parameter 

i 


Time-out = 
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Hen.cc, 


P 



1 = 




(D.4) 


We now cal culate, p 2> the probability of greater than K 
packets wna ting for Ack. If ^2,1 the traffic intensity 
at stage 2 , by analysis of a M/M/l/W system we have, 


W 

P ? = E 

j=k+i 




) p 2 ° 

W+1 j 


( 1 - 0 2 W+1 ) 


p Kt 1 / , p W-Ek 

2 u ~ 2 l 

(1- p o' ) 


U>.5) 


from (1). 2 ), (D. 4 ), (b. 5 ) we have 


(1- 


IS 


W -1 


P. 


'2 nK n K+1 




2 ( f rp 

1 1 

k=o 1 


(1- P 


W-K 


(D.6) 


For W = oo , a simple closed from solution is 
possible . 


Tho probability of greater Shon TL packets waiting, 
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or, 


P~ = 2 ( 1 - PJ pJ 

d 3 =k+ 1 




(D.8) 




^ p- Pisj-n-qHD- <§uc. 



